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Abstract 

t 

The  flow  over  an  NACA  0015  airfoil  undergoing  a  constant 
rate  of  change  of  angle  of  attack  was  experimentally  studied 
over  a  range  of  tunnel  speeds  and  rotation  rates.  Surface- 
pressure  transducers  coupled  with  a  microcomputer-based  data 
acquisition  system  were  used  to  collect  surface-pressure  data 
at  the  rate  of  4000  samples  per  second,  /data  reduction  was 

also  microcomputer-based.  SEhe'data  was  reduced  in  two  forms,  i 

jL  "...A'  'j  !••••'  •- 

— P±ts-fc>  versus  a  curves  through  stall  were  determined  for 

each  dynamic  experimental  configuration.  This  was  accom¬ 
plished  by  numerical  integration  of  the  pressure  data  at  a 
number  of  angles  through  stall,  each  data  point  representing 

the  average  of  five  experiments  at  the  same  experimental  con- 

I  o  '  '  1  j 

ditions.  These  curves  indicated  a  slight  decrease  in~C  -  a 

<j>.v  \  ■  f  >  *•  £  -j  •>  i-  .  ( "  )  ' 

slope  with  increasing  X.  b  Sec ond-lyy»%h^>  increase  in  stall 
angle  of  attack  of  the  dynamic  over  the  static  case  was  feLso 
plotted  against  a  nondimens ional  angular  rate  parameter  (de¬ 
fined  as  the  product  of  one-half  the  chord  length  and  angular 
rotation  rate,  divided  by  the  freestream  velocity) .  This  com¬ 
parison  gave  rise  to  an  apparently  universal  curve  of  non- 
dimensional  angular  rotation  rate  versus  increase  in  stall 
angle  of  attack.  This  curve  was  in  agreement  in  some  sense 
with  previous  experiments  using  stall  indicators  other  than 

the  actual  stall,  v  Data  was  collected  in  the  range  of  non- 

\ 

dimensional  angular  rates  between  .006  and  .032. 
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CONTINUED  EXPERIMENTAL  INVESTIGATION  OF  DYNAMIC  STALL 


I.  Introduction 


Background 

Dynamic  stall  is  a  physical  phenomenon  that  occurs  when 
an  airfoil  undergoes  an  uninterrupted,  dynamic  rotation 
through  its  static-stall  angle  of  attack.  Over  a  large 
range  of  angle-of-attack  rate  and  freestream  velocity  combi¬ 
nations,  the  lift  curve  continues  to  increase  beyond  the 
static-stall  point.  Though  this  is  a  transient  event,  the 
momentary  increase  in  maximum  unstalled  angle  of  attack 
yields  a  corresponding  increase  in  the  lift  generated  by  the 
airfoil.  This  greater  lift  is  of  sufficient  magnitude  to 
render  the  dynamic  stall  effect  of  some  possible  practical 
use,  and  therefore  worthy  of  further  investigation. 

The  first  formal  investigation  of  dynamic  stall  was  per¬ 
formed  by  Max  Kramer  in  1932,  and  was  prompted  by  pilot  re¬ 
ports  of  unexplained  high  lift  values  occurring  in  turbulent 
air  (Ref  1:1).  In  his  experiment,  Kramer  mounted  a  wing  on 
a  balance  in  the  test  section  of  a  wind  tunnel.  He  then 
used  movable  guide  vanes,  located  upstream  of  the  wing,  to 
produce  a  rotating  freestream  flow  in  the  test  section. 

The  resulting  angle  of  attack  experienced  by  the  wing  ranged 
from  0  to  30  degrees  (Ref  1:2-3).  Kramer  repeated  the  ex¬ 
periment  for  three  wings:  the  first  two  were  Gottingen  459 
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airfoil  cross-sections  (a  symmetric  airfoil) ,  with  different 
chord  lengths,  and  the  third  was  a  Gottingen  398  airfoil 
cross-section  (a  cambered  airfoil) . 

Kramer's  results  showed  a  direct  relationship  between 
maximum  lift  coefficient  and  the  rate  of  angle-of-attack 
change, a  ,  and  an  inverse  relationship  to  the  test  section 
velocity,  Vo.  He  introduced  a  non-dimensional  angular  rate 

Cr# 

parameter  — ,  where  c  is  the  chord  length  of  the  airfoil. 

CO 

This  parameter  allowed  Kramer  to  collapse  all  his  data  onto 
a  single  curve  given  by: 

C,  =  C1  +  0.36  (1) 

max  dyn  max  st  <*> 


Though  various  attempts  to  model  Kramer’s  experiment 
have  been  made,  the  first  successful  model  was  developed  in 
1980  by  R.  G.  Docken,  Jr.,  E.  J.  Jumper  and  J.  E.  Hitchcock 
(Ref  2) .  Their  method  involved  applying  the  momentum- integral 
equation  to  a  transient  boundary  layer  control  volume  in  un¬ 
steady  potential  flow.  The  method  was  then  applied  to  the 
case  of  an  11%  thick  symmetrical  Joukowski  airfoil  (J011) , 
and  yielded  the  result,  for  a  in  radians  per  second: 


“stall  dyn  ’  “stall  st  +  0-096  Hr* 


(2) 


As  it  stands,  Eq.  2  is  not  directly  comparable  to  Eq.  1.  It 
can  be  transformed  to  an  equivalent  version,  which  is  di¬ 
rectly  comparable,  by  assuming  that  the  only  difference  be¬ 
tween  the  static  and  dynamic  lift  curves  is  the  increased 
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stall  angle  of  attack  for  the  dynamic  case.  This  implies 
that  both  static  and  dynamic  lift  curves  have  the  same  slope. 
Correcting  this  slope  for  the  aspect  ratio  of  Kramer's  wing, 
Eq.  3  is  obtained,  which  is  directly  comparable  to  Eq.  1 
(Ref  2:3) : 

C,  =  C  +  0.301  ~  (3) 

max  dyn  max  st  « 

Comparison  of  Eq.  1  and  Eq.  3  reveals  a  remarkable 
agreement  between  Kramer's  experimental  findings  and  the 
analytical  result  of  Docken,  et.  al.  The  fact  that  Kramer 
used  a  Gottingen  459  and  Docken  a  J011  detracts  little,  if 
at  all,  from  the  credibility  of  the  comparison,  since  both 
are  symmetrical  airfoils  and  have  similar  leading  edge  geom¬ 
etries  . 

Since  the  Kramer  experiment,  which  was  the  original 
investigation  in  dynamic  stall,  a  great  deal  of  research, 
both  analytical  and  experimental,  has  been  done  in  this  area. 
However,  the  preponderance  of  work,  unlike  Kramer's,  has 
been  for  the  case  of  an  airfoil  undergoing  a  dynamic  angle- 
of-attack  change  in  a  constant-direction  freestream,  and 
most  of  this  work  has  been  done  for  the  case  of  a  sinusoi¬ 
dally  oscillating  airfoil.  The  reasons  for  this  are  fairly 
obvious.  Situations  where  dynamic  stall  research  would 
likely  find  an  application  are  those  such  as  helicopter  blades, 
turbomachinery,  and  aircraft  wing-flutter,  and  the  angle  of 
attack  variations  displayed  in  cases  such  as  these  are  likely 
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to  be  sinusoidal,  or  at  least  approximately  so. 

The  advent  of  digital  flight  control  systems, 
however,  promises  an  application  for  the  case  of  an  airfoil 
which  undergoes  an  angle  of  attack  variation  that  can  be 
described  by  a  ramp  function.  A  distinct  advantage  of  the 
ramp  angle-of-attack  variation  is  the  comparative  ease  and 
physical  clarity  with  which  a  mathematical  model  may  be 
developed,  as  opposed  to  the  case  of  the  sinusoidally  oscil¬ 
lating  airfoil.  The  mathematical  model  for  the  sinusoidal 
case  lies  in  the  realm  of  full  Navier-Stokes  solutions,  and 
amounts  to  little  more  than  a  numerical  experiment.  While 
such  an  approach  succeeds  fairly  well  in  modelling  the  re¬ 
sults  of  a  corresponding  experiment,  the  sheer  mathematical 
complexity  overwhelms  any  attempt  to  truly  understand  the 
physics  of  the  phenomenon. 

In  1979,  Deekens  and  Kuebler  (Ref  3)  undertook  an  in¬ 
vestigation  of  dynamic  stall  in  which  the  effects  of  constant 
airfoil  angular  rate  on  dynamic-stall  angle  of  attack  were 
evaluated.  Smoke-trace  flow  visualization  was  employed  in 
conjunction  with  simultaneous  movie  filming  in  order  to 
characterize  the  dynamic  stall  phenomenon  on  an  NACA  0015 
airfoil,  which  was  rotated  about  its  midchord  in  a  constant- 
velocity  freestream.  They  concluded  that  the  increase  in 
unstalled  angle  of  attack  above  that  for  static  stall  was 
dir^tly  related  to  airfoil  angular  rate,  and  inversely 
related  to  the  freestream  velocity.  Moreover,  they  were 
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able  to  predict  the  dynamic  stall  angle  of  attack  within  the 
scope  of  their  experiment,  which  was  the  Reynolds  number 
range  between  14,500  and  32,500. 

Introducing  the  nondimens ional  angular  rate  parameter, 
as  did  both  Kramer  and  Docken,  et.  al.,  Deekens  and  Kuebler 
were  able  to  collapse  their  data  onto  a  single  curve  given 
by: 


“stall  dyn 


“stall  st  + 


2.5 


l/2coi 


(4) 


A  plot  of  these  results,  showing  dynamic  stall  angle  of  at¬ 
tack  as  a  function  of  nondimensional  rotation  rate  parameter 
is  shown  in  Fig.  1,  on  the  following  page. 

By  now  introducing  the  same  assumptions  and  corrections 
as  were  used  previously  in  transforming  Eq.  2  into  Eq.  3, 

Eq.  5  is  arrived  at,  which  is  directly  comparable  to  Eq.  1 
and  Eq.  3. 


Cx  =  C  +  4.8  (5) 

max  dyn  max  stat  « 

Comparing  Eq.  1  and  Eq.  3  with  Eq.  5,  it  is  immediately 
obvious  that  the  dynamic  lift  curve  slope  implied  by  Deekens 
and  Kuebler  is  significantly  greater  than  that  given  by  the 
work  of  either  Kramer  or  Docken,  et.  al.  An  error,  either 
experimental  or  computational,  of  such  a  magnitude  to  explain 
this  apparent  discrepancy  can  be  ruled  out  since  the  results 
of  Deekens  and  Kuebler  are  substantiated  by  the  work  of 
Francis  (Ref  4),  and  by  Scheubel,  as  well  (Ref  5:1-4).  In 
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DKKKRNS  AND  KURBLRR  DATA 


Summary  of  Deekens'  ard  Kuebler's  Results 


addition,  Kramer's  work  seems  to  have  been  verified  in  an 
experiment  mentioned  by  Scheubel  (Ref  5:1). 

At  this  point  it  might  be  profitable  to  pause  for  a 
moment  to  point  out  and  emphasize  an  important  distinction. 
With  Kramer’s  experiment  configured  as  previously  described, 
the  airfoil  was  fixed  in  inertial  space  and  encountered  a 
gust.  Therefore,  a  mathematical  model  of  the  flow  over  the 
airfoil  could  justly  assume  a  Newtonian,  or  nonaccelerating, 
control  volume.  However  for  a  constant-velocity  freestream 
with  the  airfoil  being  rotated  dynamically  therein,  the  air¬ 
foil  is  rotating  with  respect  to  inertial  space.  In  this 
case,  mathematical  analysis  of  the  flow  over  the  airfoil 
cannot  be  carried  out  using  a  Newtonian  control  volume.  The 
previously  mentioned  order-of -magnitude  disagreement  between 
Kramer's  results  and  those  of  Deekens  and  Kuebler  could, 
therefore,  be  due  to  the  effect  of  the  accelerating  control 
volume . 

The  results  of  Deekens  and  Kuebler  were  recently  re¬ 
confirmed  by  Daley.  Like  Deekens  and  Kuebler,  Daley  rotated 
an  NACA  0015  airfoil  section  about  its  midchord  at  a  constant 
angular  rate  in  a  constant-velocity  freestream.  He  also  used 
smoke-trace  flow  visualization  along  with  high-speed  cinemat¬ 
ography  as  a  medium  for  recording  and  later  analyzing  his 
results.  Daley  added  a  new  dimension  to  the  experiment,  how¬ 
ever,  by  embedding  four  piezo-resistive  pressure  transducers 
in  the  airfoil  around  its  quarter-chord.  This  modification 
enabled  him  to  simultaneously  gather  two  types  of  data  re- 
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garding  the  same  phenomenon.  Using  both  movies  and  elec¬ 
tronically-gathered  pressure  information,  Daley  possessed 
an  extremely  accurate  and  sensitive  indicator  of  flow  separa¬ 
tion  at  the  quarter-chord.  Adopting  quarter-chord  flow 
separation  as  the  criterion  for  stall,  he  then  proceeded  to 
verify  a  major  portion  of  the  results  obtained  by  Deekens 
and  Kuebler.  He  also  extended  the  range  of  results  into  the 
region  of  lower  nondimensional  angular  rate,  as  shown  in 
Fig.  2,  and,  at  the  same  time,  expanded  the  Reynolds  number 
range  of  the  experiment. 

Objectives 

The  previous  research  in  dynamic  stall  for  the  constant 
rate  angle-of-attack  change  provided  a  foundation  on  which 
to  base  further  investigation.  In  view  of  this  previous 
work,  the  objectives  of  this  experiment  were  as  follows: 

1.  Equip  an  NACA  0015  airfoil  with  16  piezo-resistive 
pressure  transducers,  and  integrate  this  with  an 
automated,  microcomputer-controlled  data  collection 
system  for  gathering  time,  angle-of-attack,  and 
pressure  information. 

2 .  Develop  a  FORTRAN  software  package  to  automate  both 
the  acquisition  and  reduction  of  time,  angle-of- 
attack,  and  pressure  data. 

3.  Verify  the  operation  of  this  system,  and  then  use 
it  to  conduct  a  parametric  investigation  of  the  ef¬ 
fects  of  airfoil  angular  rate  and  freestream  veloc¬ 
ity  on  dynamic  stall. 
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DATA  SUMMARY 


Figure  2.  Data  Summary  of  Daley’s  Results 


The  following  theory  and  approach  section  is  composed 

of  six  subsections.  Each  of  these  subsections ,  in  turn, 

presents  a  brief  discussion  of  the  way  in  which  existing 

dynamic  stall  theory  or  previous  research  influenced  the  ex- 
# 

perimental  approach  in  the  current  investigation.  The  first 
subsection  provides  a  more  detailed  qualitative  description 
of  dynamic  stall  than  was  given  in  the  introduction.  The 
second  discusses  the  determination  of  pressure  coefficients 
using  ensemble-averaged  data.  Then,  the  third  subsection 
covers  discretization  of  the  pressure  distribution  defined 
by  these  pressure  coefficients,  while  the  fourth  describes 
the  integration  of  this  discretized  pressure  distribution. 

The  fifth  subsection  considers  the  computation  of  force  co¬ 
efficients  using  the  results  of  the  integration,  and  the 
sixth  presents  a  brief  narrative  concerning  the  problem  of 
data  acquisition. 

Dynamic  Stall  in  Contrast  to  Static  Stall 

Stall,  whether  static  or  dynamic,  is  defined  as  having 
occurred  when  the  boundary  layer  has  separated  from  the  air¬ 
foil  to  such  a  degree  that  the  resulting  pressure  distribu¬ 
tion  ceases  to  yield  any  increase  in  lift  for  further  angle- 
of-attack  increase.  Obviously,  the  boundary  layer  inter¬ 
actions  for  static  and  dynamic  stall  must  differ  significant- 


ly  to  produce  the  dramatic  dissimilarities  previously  men¬ 
tioned  in  the  introduction. 

In  the  familiar  case  for  static  stall,  a  boundary  layer 
under  the  influence  of  an  adverse  pressure  gradient  will 
eventually  separate  from  the  surface  of  the  airfoil  at  the 
point  where  the  shear  stress  at  the  wall  vanishes.  The 
point  where  flow  separates  and  where  it  reverses  are  coinci¬ 
dent  for  static  stall.  Thus,  the  wake  formed  by  this  viscous 
interation  is  large,  and  appreciably  distorts  the  potential 
flow  field  around  the  airfoil.  For  the  static  case,  the 
stall  angle  of  attack  is  relatively  constant,  being,  at  most, 
a  weak  function  of  Reynolds  number  (Ref  6:248). 

In  dynamic  stall,  as  with  static  stall,  the  boundary 
layer  under  the  influence  of  an  adverse  pressure  gradient 
also  eventually  reaches  the  condition  of  vanishing  wall 
shear.  The  similarity  ends  here,  however,  as  the  point  of 
reversed  flow  is  no  longer  coincident  with  the  point  of  zero 
wall  shear,  but  is  delayed  for  some  distance  downstream  of 
the  zero-shear  point.  The  point  of  separation  for  the  case 
of  dynamic  stall  is  determined  by  the  MRS  (Moore-Rott-Sears) 
criterion  (Ref  7:113-144).  This  difference  substantially 
reduces  the  wake  size  and  corresponding  potential  flow  field 
distortion  in  comparison  to  the  static  case  (Ref  8:2-3).  In 
addition  to  the  MRS  separation  criterion,  other  effects  ap¬ 
pear  to  be  at  work,  as  well  (Ref  15).  Thus,  it  is  clear 
that  the  dynamic  stall  angle-of-attack  is  a  complex  function 
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of  freestream  velocity,  airfoil  angular  rate,  and  even  air¬ 
foil  section  geometry. 

The  fact  that  the  wake  size  and  corresponding  potential 
flow  field  distortion  associated  with  dynamic  stall  are 
small  relative  to  their  static-stall  counterparts  is  favor¬ 
able  to  the  present  study.  Tunnel  wall  interference  effects 
are  therefore  correspondingly  small  compared  to  those  en¬ 
countered  in  the  same  flow  regime  for  steady-state  phenomena. 
This  implies  that  the  effects  of  streamline  curvature  and 
wake  blockage  can  probably  be  considered  negligible  in  dy¬ 
namic  stall  testing.  A  more  detailed  discussion  of  these 
test  section  wall  effects  is  presented  in  Appendix  A. 

Determination  of  Pressure  Coefficients 

Because  of  inevitable  freestream  irregularities  in  the 
test  section,  pressure  measurements  taken  at  the  same  loca¬ 
tion  on  the  airfoil,  but  at  different  times,  are  not  constant. 
These  irregularities  can  be  filtered  out  using  ensemble¬ 
averaging  which,  at  the  same  time,  preserves  those  pressure 
fluctuations  due  only  to  the  dynamic  stall  phenomenon.  For 
this  experiment,  pressure  data  from  five  different  airfoil 
rotations  at  the  same  angular  rate  and  freestream  velocity 
were  averaged  to  yield  an  ensemble-averaged  data  set. 

Using  this  ensemble-averaged  pressure  data,  a  method  of 
calculating  the  pressure  coefficient  at  any  chord  location 
on  the  airfoil  must  be  found.  This  method  should  use  as  its 
inputs  physical  parameters  which  can  be  readily  sensed  or 
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measured.  The  standard  equation  for  the  computation  of 
pressure  coefficient  is  given  by: 


c 


P 


^loc  P°° 

1/2  Pa> 


(6) 


where  Ploc  is  the  local  static  pressure  at  some  location  on 
the  airfoil,  p  is  the  freestream  static  pressure,  and  p 
and  are  the  freestream  density  and  velocity,  respectively. 
The  local  pressure  anywhere  on  the  airfoil  can  be  expressed 
as : 


ploc  -  4  ptran  +  pa  17 

where  Ploc  retains  the  same  definition  it  had  in  Eq.  6,  pg 
is  some  reference  pressure,  and  <4Ptran  is  the  differential 
pressure  between  these  two.  Substituting  Ploc  from  Eq.  7 
into  Eq.  6  yields  the  relationship: 
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(A  p 


tran 


1/2p. 


Pa5  ‘ 


(8) 


Regrouping  the  terms  in  the  numerator  and  noting  that  the 
denominator  is  equivalent  to  pQ  -  pa  under  the  incompressible 
Bernoulli  equation,  Eq.  8  becomes: 


A  P 
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P  J 


Po  -  Pa 


(9) 


Eq.  9  requires  that  three  quantities  be  known  in  order 

to  determine  the  corresponding  pressure  coefficient.  The 

first,  Ap  ,  is  the  difference  in  pressure  between  some 
trsn 
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constant  reference  pressure,  p  .  and  the  pressure  at  a  cer- 

a 

i 

tain  point  on  the  surface  of  the  airfoil.  This  differential 

I 

pressure  was  sensed  by  a  transducer  mounted  in  the  airfoil.  j 

The  second,  p^  -  p^,  is  the  pressure  difference  between  the 

reference  pressure  and  test  section  static  pressure,  while 

the  third,  pQ  -  pB,  is  the  pressure  difference  between  test 

section  stagnation  pressure  and  static  pressure.  Since  the 

reference  pressure  must  be  easily  accessible  as  well  as 

constant,  ambient  room  pressure  constituted  a  good  choice, 

although  any  constant,  easily  accessible  pressure  source 

would  have  been  acceptable. 

Discretization  of  the  Pressure  Distribution 

The  mathematical  procedure  developed  in  the  preceding 
subsection  facilitates  pressure  coefficient  determination  at 
any  point  on  the  airfoil  where  a  transducer  is  located. 

Hence,  to  minimize  the  error  inherent  in  discretizing  the 
dynamic  stall  pressure  distribution,  two  basic  issues  must 
be  addressed.  First,  an  acceptable  number  of  transducers 
must  be  determined  and,  second,  the  optimum  distribution  of 
these  transducers  has  to  be  established. 

Obviously,  a  greater  number  of  transducers  yields  a  re¬ 
duced  discretization  error.  However,  an  upper  limit  on  this 
number  is  eventually  reached,  and  is  dictated  mainly  by 
practical  considerations.  In  this  experiment,  16  transducers 
were  used,  the  number  being  limited  chiefly  by  the  capacity 
of  the  electronic  data  acquisition  system  (c.f.  below). 
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Placement  of  these  16  transducers  on  the  airfoil  was 
governed  by  the  requirement  that  the  resulting  discretized 
pressure  distribution  portray  as  accurately  as  possible  the 
actual  pressure  distribution.  Thus,  the  transducers  were 
concentrated  in  the  region  of  the  airfoil  where  the  dynamic 
stall  pressure  distribution  was  anticipated  to  have  the 
largest  gradient.  The  pressure  distribution  obtained  by 
McAlister,  Carr  and  McCroskey  (Ref  9:51)  for  an  NACA  0012 
airfoil  oscillating  at  a  nondimensional  angular  rate  com¬ 
parable  to  that  for  the  present  study  was  useful  here  as  a 
guide.  Accordingly,  the  transducers  were  distributed  most 
densely  on  the  upper  surface  of  the  airfoil  near  the  leading 
edge.  Fig.  3  illustrates  the  transducer  placement  chosen  to 
minimize  the  error  inherent  in  discretizing  the  anticipated 
pressure  distribution. 

As  shown  in  Fig.  3,  there  is  no  transducer  at  the 
trailing  edge  of  the  airfoil  and,  as  a  consequence,  no  way 
to  determine  the  pressure  at  the  trailing  edge.  However, 
the  results  of  McAlister,  et.  al.,  were  obtained  using  an 
airfoil  that  had  a  pressure  transducer  located  at  the  98 
percent  chord  position  (Ref  9:51).  Again  using  their  results 
(Ref  9:51)  as  a  guide,  it  appears  that  the  trailing  edge 
pressure  coefficient  can  be  approximated  in  this  experiment 
by  using  a  linear  curve  fit  involving  the  rearmost  two  trans¬ 
ducers  on  the  upper  surface  of  the  airfoil. 
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Figure  3 •  Upper  and  Lower  Surface  Pressure 
Transducer  Locations 


Integration  of  the  Pressure  Distribution 


A  numerical  integration  was  performed  on  the  discretized 
pressure  distribution  to  obtain  the  corresponding  force  co¬ 
efficients.  McAlister,  et.  al.,  found  that  cubic  and  variable 
power  splines  applied  to  the  discrete  data  points  did  not 
yield  acceptable  accuracy.  The  spline  fits  caused  large 
overshoots  that  made  this  method  unsatisfactory  in  general 
application  (Ref  9:12).  Therefore,  following  McAlister, 
et.  al.,  all  integrations  in  this  investigation  were  based 
on  the  trapezoidal  rule. 

Determination  of  Force  Coefficients 

The  body  axis  normal  force  can  be  decomposed  into  its 
constituent  lift  and  drag  forces  in  the  wind  axis  system, 
and  the  drag  force  can  be  further  separated  into  pressure 
drag  and  skin-friction  drag.  The  integration  of  the  pressure 
distribution  properly  accounts  for  the  lift  and  pressure  drag 
forces,  but  is  incapable  of  predicting  the  skin-friction  drag 
(a  viscous  effect).  This  is  not  a  serious  source  of  error, 
however,  for  the  following  reasons. 

At  low  angles  of  attack,  where  flow  remains  attached 
over  most  of  the  airfoil  and  skin  friction  is  the  major 
source  of  drag,  the  angle  between  the  normal  force  vector 
and  the  drag  force  vector  is  large.  Since  this  angle  is 
large,  the  projection  of  the  drag  force  vector  onto  that  of 
the  normal  force  is  correspondingly  small,  thereby  making 
any  error  in  the  drag  force  negligible  with  respect  to  the 
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magnitude  of  the  normal  force  vector. 


At  higher  angles  of  attack,  the  projection  of  the  drag 
force  vector  onto  the  normal  force  vector  and  the  normal 
force  vector,  itself,  are  of  comparable  magnitude.  In  this 
case,  however,  the  flow  is  no  longer  fully  attached  to  the 
airfoil.  As  a  result,  skin-friction  drag  is  overshadowed  by 
pressure  drag.  Thus,  only  the  drag  force  coefficient  would 
be  significantly  in  error,  and  even  then,  only  at  low  angles 
of  attack. 

To  accurately  take  into  account  skin-friction  drag, 
direct  force  measurement  could  be  used.  The  preferred 
method,  however,  still  consists  of  integrating  the  pressure 
distribution  on  the  airfoil.  Proceeding  in  this  manner  not 
only  yields  the  normal  force  coefficient,  but  characterizes 
the  pressure  distribution  as  a  necessary  intermediate  step. 
This  characterization,  in  turn,  constitutes  an  indispensable 
tool  in  the  analysis  of  the  flow  over  the  airfoil  as  it 
undergoes  dynamic  stall. 

The  Problem  of  Data  Acquisition 

Measurement  of  the  physical  parameters  associated  with 
dynamic  stall  presents  a  unique  problem,  due  to  the  transient 
nature  of  the  phenomenon.  The  measurement  system  used  must 
be  not  only  accurate,  but  fast,  as  well.  The  solution  to 
this  problem  has  taken  many  forms,  beginning  with  Kramer, 
who  used  a  balance  system  to  measure  the  aerodynamic  forces 
on  the  wing  as  the  freestream  flow  was  rotated  past  it. 
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Deekens  and  Kuebler  used  high-speed  cinemaphotography  of 
smoke  traces  to  ascertain  airfoil  rotation  rate  and  the 
angle  of  attack  at  which  dynamic  stall  occurred.  Daley, 
again,  used  movies  of  smoke  traces,  but  simultaneously 
gathered  digital  electronic  position  and  pressure  data  using 
four  electronic  transducers  embedded  in  the  quarter-chord 
region  of  the  airfoil.  McAlister,  et.  al.,  used  an  airfoil 
equipped  with  16  pressure  transducers,  and  collected  analog 
electronic  position  and  pressure  data. 

In  the  present  study,  digital  position  and  pressure  in¬ 
formation  was  collected  using  an  airfoil  instrumented  with 
16  pressure  transducers  (c.f.  above).  To  gather  digital  in¬ 
formation  for  the  purpose  of  later  constructing  dynamic  lift 
curves,  it  was  insured  that  the  sample  rate  was  high  enough 
to  yield  the  resolution  required  for  such  a  purpose.  The 
lower  threshold  on  sample  rate,  in  this  case,  was  approxi¬ 
mately  300  data  samples  per  second  (Ref  10:7).  The  data  ac¬ 
quisition  system  used  in  the  present  study  had  the  capability 
to  meet  and  exceed  this  criterion  by  a  wide  margin  (c.f. 
below) . 
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III.  Facilities  and  Instrumentation 


Smoke  Tunnel 

The  present  study  was  conducted  in  the  AFIT  smoke 
tunnel  located  in  Building  640,  Area  B,  Wright-Patterson 
AFB ,  Ohio.  The  test  section  is  59  inches  long,  39.5  inches 
high,  and  2.75  inches  deep.  This  facility  is  capable  of 
test  section  velocities  between,  approximately,  10  and  50 
feet  per  second.  The  smoke  tunnel  and  its  capabilities  are 
further  described  by  Sisson  (Ref  11) ,  and  Baldner  (Ref  12) . 

Velocity  Measurement 

Test  section  static  and  total  pressure  were  measured 
using  a  standard  hemispherical-head  Pitot-static  probe  in 
conjunction  with  a  Meriaro  A-937  water  micromanometer.  These 
pressures  were  used  to  establish  the  test  section  velocity 
during  data  collection  and,  later,  to  determine  the  pressure 
coefficients  when  reducing  the  data.  This  made  the  accurate 
measurement  of  these  pressures  a  crucial  step  in  the  experi¬ 
ment.  It  was  therefore  necessary  to  locate  the  probe  as 
close  to  the  airfoil  as  possible,  without  creating  mutual 
interference  between  the  two,  due  to  the  presence  of  a  sig¬ 
nificant  longitudinal  static  pressure  gradient  (see  Appendix 
B)  . 

Considering  this,  the  probe  hole  was  drilled  into  the 
floor  of  the  tunnel  at  a  location  31  inches  downstream  of 
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the  point  where  the  test  section  begins.  With  the  airfoil 
at  the  zero  angle-of-attack  position,  the  tip  of  the  probe 
was  directly  below  the  leading  edge  of  the  airfoil,  and  the 
probe  static  ports  were  coplanar  with  the  13  percent  chord 
station  on  the  airfoil. 

Airfoil 

The  NACA  0015  airfoil  used  in  this  experiment  measured 
12.2  inches  chord  and  2.63  inches  span.  It  consisted  of  a 
hollow  mahogany  shell  closed  on  both  sides  by  aluminum  end- 
plates,  which  were  sealed  to  the  shell  with  silicone  rubber 
adhesive  sealant.  To  reduce  the  frictional  effects  of  the 
tunnel  walls  on  the  rotating  airfoil.  Teflon  bearing  sur¬ 
faces  were  fitted  to  both  the  front  and  rear  endplates.  The 
rear  endplate  was  rigidly  attached  to  a  14  inch  tubular 
aluminum  shaft  with  an  outside  diameter  of  .75  inches.  This 
aluminum  shaft  was  slotted  at  its  midpoint  to  admit  ambient 
atmospheric  pressure  to  the  interior  of  the  airfoil,  as  well 
as  the  transducer  leads.  Sixteen  ports  were  bored  into  the 
wooden  shell  using  a  number  38  drill,  to  accept  the  16 
transducers  (c.f.  below). 

Transducers 

The  transducers  chosen  for  use  in  this  experiment  were 
Endevco  8506-2  and  8507-2  miniature  piezo-resistive  pressure 
transducers,  the  only  difference  between  the  8506  and  8507 
being  that  the  8506  had  a  threaded  mounting.  Both  types  of 
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transducers  had  a  maximum  range  of  plus  or  minus  two  psig, 
and  required  an  excitation  voltage  of  10.00  volts  DC.  Ex¬ 
citation  voltage  was  provided  by  a  Hewlett-Packard  6205B 
Dual  DC  Power  Supply,  and  monitored  by  a  Hewlett-Packard 
34701A  DC  voltmeter  with  a  34740A  digital  display  insert, 
which  allowed  voltage  to  be  read  to  three  decimal  places. 
Resonance  frequency  for  both  types  of  transducers  was  45,000 
Hertz.  Thus,  the  frequency  response  of  the  transducer  had 
a  negligible  effect  on  the  results  obtained  in  this  experi¬ 
ment  . 

The  transducers  were  mounted  in  the  ports  in  the  airfoil 
according  to  the  specifications  provided  by  Endevco  (Ref  13) . 
General  Electric  RTV  silicone  rubber  adhesive  sealant  was 
used  here  as  a  bonding  agent.  After  completing  the  electri¬ 
cal  connections  between  the  transducers  and  the  microcomputer, 
the  transducers  were  recalibrated  (c.f.  below).  In  many 
cases,  the  sensitivities  obtained  in  the  recalibration  were 
significantly  different  from  those  provided  by  Endevco. 

Drive  Mechanism 

The  airfoil  was  rotated  using  a  TRW  Globe  Model  5A2298-4 
12  volt  DC,  constant-speed  planetary  gearmotor  with  a  525:1 
reduction  ratio.  The  attributes  of  this  unit  gave  it  two 
definite  advantages  over  the  alternatives  considered.  First, 
its  power  requirements  could  be  met  using  a  standard  labora¬ 
tory  bench  power  supply,  which,  in  this  experiment,  was  the 
Hewlett-Packard  6205C  Dual  DC  Power  Supply.  By  adjusting 
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the  supply  voltage  between  6  and  12  volts,  the  rotation  rate 
at  the  output  shaft  could  be  varied  between  approximately  30 
and  95  degrees  per  second.  Second,  the  high  reduction  ratio 
gave  a  high  output  torque.  This  high  output  torque,  in  turn, 
spun  the  output  shaft  up  to  constant  speed  in  less  than  .01 
seconds,  which  was  negligible  compared  to  the  time  required 
to  rotate  the  airfoil  through  the  dynamic  stall  angle  of 
attack.  Thus,  airfoil  angular  rate  was  easily  governed  and, 
once  the  motor  voltage  was  set,  the  variation  in  the  output 
angular  rate  was  negligible.  The  test  case  was  therefore 
easily  repeatable.  A  spring-loaded  double-pole,  double-throw 
toggle  switch  was  used  to  rotate  the  airfoil,  stop  it,  and 
then  return  it  to  the  initial  position  in  preparation  for 
the  next  run. 

The  airfoil  angle-of-attack  transducer  consisted  of  a 
Helipot  #7216  1000  ohm,  ten-turn  potentiometer.  The  potentio¬ 
meter  was  coupled  to  the  airfoil  shaft  through  a  gear  train 
having  a  33:1  ratio,  giving  a  full  ten  turns  at  the  potentio¬ 
meter  for  109.1  degrees  airfoil  rotation.  The  potentiometer 
was  excited  at  10  volts,  simultaneously  using  the  same  volt¬ 
age  source  as  used  for  transducer  excitation,  and  its  output 
was  fed  to  one  of  two  analog-to-digital  conversion  boards  in 
the  computer  (c.f.  below). 

Data  Acquisition  System 

The  microcomputer  system  consisted  of  a  Heathkit  H-19 
terminal,  a  Tarbell  Model  VDS-IID  dual  eight-inch  floppy  disk 
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drive,  and  an  Electronic  Control  Technology  S-10G  bus  equipped 
with  an  SD  Systems  SBU-100  Single  Board  Computer,  an  SD  Sys¬ 
tems  Expandoram  II  board,  and  an  MD2022  Tarbell  Disk  Control¬ 
ler  board.  To  perform  its  digital  data  gathering  function, 
this  system  was  augmented  with  two  Dual  Systems  Control 
Corporation  AIM-12  analog  input  module  boards. 

The  AIM-12  is  a  high  speed,  multiplexed  analog-to-digital 
data  acquisition  module  compatible  with  the  standard  S-100 
bus.  The  analog-to-digital  conversion  subsystem  on  the  board 
can  be  operated  in  one  of  two  modes.  The  unipolar  mode  re¬ 
quires  that  the  input  voltage  to  the  A/D  converter  be  between 
0  and  10  volts,  while  the  bipolar  mode  accepts  input  voltages 
from  -5  volts  to  +5  volts. 

Since  the  analog  signal  to  be  digitized  will  usually 
not  give  full-scale  deflection  in  either  of  these  modes,  the 
board  has  a  preconditioning  subsystem  consisting  of  a  multi¬ 
plexed,  precision  instrumentation  amplifier.  By  changing 
the  resistance  in  the  feedback  loop  of  the  amplifier,  gains 
from  1  to  100  can  be  selected.  Single-ended  amplifier  opera¬ 
tion  allows  32  separate  analog  inputs  to  the  multiplexer, 
while  differential  operation  limits  the  multiplexer  to  16 
inputs.  Differential  operation,  however,  takes  advantage  of 
the  high  common  mode  rejection  ratio  of  the  amplifier,  which 
is  80  decibels  at  a  gain  of  1  and  increases  to  114  decibels 
for  the  maximum  gain  of  100. 

As  mentioned  previously,  two  of  these  AIM-12  boards  were 
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used  in  the  data  acquisition  system.  The  board  responsible 
for  the  collection  and  digitization  of  the  electrical  signals 
from  the  16  pressure  transducers  was  configured  for  gain  100 
amplification  in  the  differential  mode,  and  bipolar  A/D  con¬ 
version.  The  electrical  signals  originating  at  the  pressure 
transducers  had  a  magnitude  of,  approximately,  15  millivolts. 
Although  the  gain  of  100  resulted  in  no  more  than  30  percent 
of  full-scale  on  the  A/D  converter,  the  high  common  mode  re¬ 
jection  ratio  was  very  effective  in  cancelling  noise  in  the 
system  and  the  resulting  overall  system  accuracy  was  quite 
good. 

Another  point  that  should  be  discussed  is  the  reason 
for  choosing  bipolar  A/D  conversion.  The  highest  pressure 
that  will  be  experienced  in  any  flow  condition  is  that  at 
stagnation.  Since  the  smoke  tunnel  is  of  the  drawdown  vari¬ 
ety,  the  stagnation  pressure  in  the  test  section  is  the  same 
as  ambient  pressure  if  inlet  losses  are  assumed  negligible, 
and  lower  if  they  are  taken  into  account.  Therefore,  it  is 
reasonable  to  assume  that  the  highest  pressure  that  could 
possibly  occur  in.  the  test  section  is  somewhat  lower  than 
ambient  pressure.  The  transducers  in  the  present  configura¬ 
tion  sense  the  difference  between  some  local  pressure  on  the 
airfoil  and  ambient,  or  room,  pressure.  With  the  highest 
pressure  in  the  test  section  being  somewhat  lower  than  am¬ 
bient  pressure,  it  seems  that  the  transducers  would  output 
only  negative  voltages  and  the  unipolar  A/D  conversion  mode 
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could  be  used,  thereby  cutting  the  worst  case  error  in  half. 
However,  due  to  zero  drift,  the  transducers  did,  in  some 
instances,  output  a  voltage  in  the  positive  range,  obviating 
the  need  to  use  the  bipolar  A/D  conversion  mode. 

The  other  AIM-12  board  was  responsible  for  the  collection 
and  digitization  of  the  signal  from  the  position  potentiome¬ 
ter,  which  varied  between  0  and  10  volts.  This  board  was 
configured  for  gain  1  amplification  and  unipolar  A/D  conver¬ 
sion. 
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IV .  Experimental  Procedure 

Transducer  Calibration 

All  16  transducers  in  the  airfoil  were  calibrated  prior 
to  the  first  data  collection  run.  This  calibration  procedure 
was  repeated  every  three  weeks  until  data  gathering  was  com¬ 
plete.  A  more  complete  description  of  transducer  calibration 
is  presented  in  Appendix  C. 

Data  Collection 

To  prepare  the  system  for  a  data  collection  run,  all 
three  voltmeters,  both  power  supplies,  and  the  computer  were 
allowed  to  warm  up  at  for  a  minimum  of  four  hours  before  any 
data  was  taken.  This  was  done  to  allow  any  large  electrical 
transients  in  the  system  to  die  out,  insuring  steady-state 
operation  during  data  collection. 

The  first  step  in  making  a  data  collection  run  was  to 
execute  the  data  acquisition  program,  BIG  (see  Appendix  D) . 
The  remainder  of  the  experimental  procedure  was  then  auto¬ 
matically  indicated  at  the  computer  terminal  by  requests  for 
input,  or  as  simple  instructions  concerning  equipment  opera¬ 
tion.  The  following  discussion  constitutes  a  summary  of  the 
data  collection  sequence. 

The  first  set  of  inputs  to  the  computer  consisted  of 
date,  time,  temperature,  and  barometer.  These  were  then 
echoed  back  to  the  operator  for  verification  and,  if  desired, 
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corrections  to  these  inputs  could  then  be  made  before  they 
were  written  to  the  disk  file.  Next,  the  zero-input  readings 
for  the  16  transducers  were  taken,  displayed  on  the  terminal 
screen,  and  written  to  disk.  At  this  point  the  operator  was 
instructed  by  the  computer  to  turn  on  the  tunnel  motors  and 
bring  the  test  section  flow  up  to  the  desired  velocity. 

The  next  set  of  inputs  consisted  of  the  two  different 
micromanometer  readings,  the  voltages  corresponding  to  the 
90  and  0  degree  angle-of-attack  positions,  and  the  voltage 
that  was  to  be  supplied  to  airfoil  drive  motor.  The  first 
micromanometer  reading  was  the  difference  between  ambient 
pressure  and  test  section  static  pressure.  This  was  obtained 
by  connecting  the  tube  from  the  Pitot-static  probe  static 
ports  to  one  leg  of  the  micromanometer  and  leaving  the  other 
leg  open  to  the  ambient  air.  The  second  micromanometer 
reading  was  the  difference  between  the  test  section  total 
and  static  pressures,  and  was  obtained  by  connecting  tht-  tube 
from  the  probe  total  pressure  port  to  the  other  leg  of  the 
micromanometer . 

The  voltages  corresponding  to  90  and  0  degrees  angle-of- 
attack  were  measured  using  a  digital  voltmeter  connected  to 
the  position  potentiometer  during  this  phase  of  the  data 
collection  run.  The  90  and  0  degree  angle-of-attack  positions 
were  indicated  by  markers  attached  to  the  front  glass  wall  of 
the  test  section.  After  these  voltages  were  input  to  the 
computer,  the  motor  voltage  was  input,  and  all  values  input 
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in  this  set  were  echoed  at  the  terminal  screen  for  verifica¬ 
tion.  The  operator  then  had  the  option  to  write  these  values 
to  disk  or  to  enter  them  again. 

The  next  part  of  the  program  did  the  actual  dynamic  stall 
data  collection  for  five  consecutive  airfoil  rotations  at  one 
motor  voltage  and  one  test  section  velocity.  The  following 
description  is  representative  of  the  sequence  that  the  opera¬ 
tor  performed  for  one  such  airfoil  rotation. 

The  operator  would  first  input  the  number  of  samples  to 
be  taken,  as  well  as  his  choice  of  manual  or  automatic  trig¬ 
ger.  The  number  of  samples  and  choice  of  trigger  remained 
consistent  over  the  five  consecutive  rotations  to  avoid  dif¬ 
ficulty  in  data  reduction.  After  the  airfoil  had  been  rotated 
through  dynamic  stall  and  returned  to  zero  angle  of  attack, 
the  computer  would  output,  to  the  terminal,  the  number  of 
samples  actually  taken  and  the  angular  rotation  rate  of  the 
airfoil  in  degrees  per  second.  At  this  point,  the  operator 
had  the  option  to  write  the  data  set  just  collected  to  disk, 
or  to  repeat  the  rotation  to  get  a  new  data  set.  The  rota¬ 
tion  was  repeated  if  the  rotation  rate  was  not  within  two 
degrees  per  second  of  the  angular  rates  for  the  rotations 
previously  completed  at  the  same  test  section  velocity  and 
motor  voltage.  If  it  did  meet  this  criterion,  the  data  set 
just  collected  was  written  to  disk. 

After  five  satisfactory  dynamic  stall  data  files  had 
been  obtained,  a  static  lift  curve  was  determined  for  the 
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same  test  section  velocity.  The  next  part  of  the  program  did 
this  by  first  instructing  the  operator  to  rotate  the  airfoil 
to  the  desired  static  angle  of  attack.  Then,  at  the  command 
of  the  operator,  the  transducers  were  sampled  the  specified 
number  of  times,  and  the  resulting  normal  force  coefficient 
was  computed  and  displayed  at  the  terminal .  This  value  was 
recorded  and  the  procedure  was  repeated  a  sufficient  number 
of  times  at  successively  higher  angles  of  attack  to  define  a 
static  lift  curve.  The  static  angle  of  attack  was  read  from 
a  clear  plastic  protractor  taped  to  the  front  wall  of  the 
test  section.  When  the  static  lift  curve  had  been  determined 
to  the  satisfaction  of  the  operator,  the  data  collection  pro¬ 
gram,  BIG,  was  terminated,  and  the  tunnel  shut  down  until  the 
next  run. 

Velocities  and  Reynolds  Numbers 

Using  the  procedure  outlined  above,  test  runs  were  con¬ 
ducted  at  five  test  section  velocities  between  26  and  48  feet 
per  second,  inclusive.  Although  the  smoke  tunnel  was  capable 
of  test  section  velocities  as  low  as  10  feet  per  second,  any 
data  gathered  at  velocities  below  approximately  25  feet  per 
second  was  not  of  acceptable  quality  for  two  reasons.  First, 
the  magnitude  of  the  resulting  signal  was  of  low  enough  mag¬ 
nitude  to  render  the  resulting  signal-to-noise  ratio  unac¬ 
ceptable.  Second,  the  percent  full-scale  deflection  at  the 
analog-to-digital  converter  was  small  enough  to  make  the  re¬ 
sulting  analog-to-digital  resolution  unacceptable.  At  each 
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test  section  velocity,  a  run  was  accomplished  for  each  of 
four  motor  voltages,  giving  a  total  of  20  test  runs.  The 
resulting  Reynolds  numbers,  based  on  airfoil  chordlength, 
ranged  from  1.58  x  10^  to  2.81  x  10^.  As  such,  all  data  was 
collected  in  a  flow  regime  generally  considered  laminar, 
based  on  Reynolds  number. 
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V.  Data  Reduction  and  Discussion  of  Results 


Data  Reduction 

The  raw  data  files  generated  during  the  test  runs  and 
stored  on  disk  were  later  reduced  using  the  data  reduction 
program  RED,  a  copy  of  which  is  included  in  Appendix  D. 

This  program  first  ensemble  averaged  the  five  data  files 
(all  five  of  which  were  gathered  at  one  test  section  velo¬ 
city  and  one  motor  voltage) ,  then  converted  the  averaged 
digital  voltages  to  angles  of  attack  and  pressures.  These 
pressures  were  converted  to  pressure  coefficients  which  were, 
in  turn,  used  to  construct  a  pressure  distribution.  This 
pressure  distribution  was  then  integrated  using  the  trap¬ 
ezoidal  rule,  which  gave  a  dynamic  normal  force  coefficient 
for  a  corresponding  angle  of  attack.  RED  reduced  every  tenth 
data  set,  a  data  set  being  the  elapsed  time,  position  vol¬ 
tage,  and  the  voltages  from  the  16  transducers  which  were 
collected  in  one  sampling  pass. 

Since  RED  reduced  every  tenth  data  set,  another  pro¬ 
gram  was  needed  to  better  resolve  the  point  at  which  dynamic 
stall  actually  occurred.  This  program,  DEN,  reduced  every 
data  set  within  a  specified  interval  on  either  side  of  the 
anticipated  dynamic  stall  angle  of  attack.  The  anticipated 
dynamic  stall  angle  of  attack  was  obtained  by  first  running 
RED  and  finding  the  approximate  angle  of  attack  at  which 
dynamic  stall  occurred.  A  copy  of  DEN  is  included  in  Appen- 
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dix  D.  Programs  RED  and  DEN  both  include  comments  to  facili¬ 
tate  understanding  of  the  data  reduction  routines  used. 

Discussion  of  Results 

The  reduced  data  obtained  using  RED  and  DEN  was  plotted, 
giving  a  dynamic  lift  curve  (.that  is,  for  a  f  0)  for  each 
of  the  twenty  different  combinations  of  motor  voltage  and 
test  section  velocity.  A  representative  set  of  these  plots 
for  34.4  feet  per  second  test  section  velocity  and  the  full 
range  of  four  motor  voltages  is  presented  in  Figs.  4  through 
7.  A  static  lift  curve  for  the  same  test  section  velocity 
is  also  included  in  this  set,  in  Fig.  8.  Appendix  E  con¬ 
tains  the  remainder  of  the  plotted  results  and  the  associ¬ 
ated  static  lift  curves. 

Comparison  of  the  static  lift  curve  in  Fig.  8  with  the 
one  given  in  Ref  6  for  a  NACA  0015  at  a  similar  Reynolds 
number  leads  to  the  following  three  observations.  First, 
the  slope  of  the  curve  obtained  in  this  experiment  is  slight¬ 
ly  lower  than  that  shown  in  Ref  6.  Second,  the  character 
of  the  curve  in  the  range  between  eight  and  fourteen  degrees 
angle  of  attack  differs  somewhat  from  the  one  in  Ref  6. 
Finally,  the  static-stall  angle  of  attack  is  approximately 
six  degrees  higher  than  would  be  expected  based  on  Ref  6. 

The  difference  in  lift  curve  slope  could  be  due,  in 
part,  to  the  influence  of  three-dimensional  effects  on  the 
airfoil.  One  source  of  three-dimensional  effects  was  the 


ANGLE  OF  ATTACK  (DEGREES) 


Figure  4.  Force  Coefficient  vs.  Angle  of  Attack 

for  0=32.6  Deg/Sec  and  A  Ft/Sec 

(  dND=0.009  ) 


Figure  5.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =52.7  Deg/Sec  and  ^  =  34.4  Ft/Sec 
(  ®NDs0'0l4  ) 
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ANGLE  OF  ATTACK  (DEGREES) 


Figure  6.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =74 .4  Deg/Sec  and  V  =3^.^  Ft/Sec 
(  ®nds0*°19  ) 


FORCE  COEFFICIENT 


Figure  7.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =95  .9  Deg/Sec  and  Vc_=34.4  Ft/Sec 
(aND=0.025  ) 
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Figure  8.  Force  Coefficient  vs.  Angle  of  Attack 
for  6=0.0  Deg/Sec  and  Ft/Sec 
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interaction  between  the  airfoil  and  the  boundary  layer  along 
the  walls  of  the  test  section  (see  Appendix  B) .  Another 
source  of  these  effects  was  the  clearance  between  the  air¬ 
foil  and  the  walls  of  the  test  section  necessary  to  allow 
airfoil  rotation.  This  clearance  varied  from  .125  inches 
with  the  tunnel  shut  down,  to  less  than  .05  inches  for  a 
test  section  velocity  of  48  feet  per  second.  To  minimize 
the  leakage  through  this  space,  felt  was  glued  to  both  sides 
of  the  airfoil. 

The  slight  differences  in  the  character  of  the  static 
lift  curve  at  higher  angles  of  attack  might  suggest  a  possi¬ 
ble  leading  edge  separation  bubble  (Ref  14:2).  However, 
this  difference  is  so  slight  that  it  could  be  attributed 
solely  to  data  averaging  effects  (c.f.  above).  Further, 
smoke-on  films  made  by  Daley  (Ref  10)  did  not  appear  to  in¬ 
dicate  the  presence  of  such  a  bubble  up  to  stall,  at  least 
for  the  dynamic  cases.  Since  Daley  filmed  only  the  stall 
location  for  the  static  case,  it  can  only  be  stated  that  a 
separation  bubble  was  apparently  not  present  in  the  dynamic 
cases. 

Another  factor  possibly  contributing  to  the  difference 
in  the  static  lift  curve  is  tunnel  blockage  due  to  the  air¬ 
foil  itself  as  well  as  the  wake  it  generates  at  the  higher 
angles  of  attack.  An  explanation  of  these  blockage  effects 
is  included  in  Appendix  A. 
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Finally,  it  is  difficult  to  formulate  a  hypothesis  as 
to  why  the  stall  angle  of  attack  in  the  present  study  was 
higher  than  that  given  for  the  same  airfoil  in  Ref  6.  High¬ 
er  effective  Reynolds  number,  whether  due  to  turbulence  in 
the  test  section  or  airfoil  roughness,  can  probably  be  ruled 
out,  though.  While  there  was  an  increase  in  the  angle  of 
attack  at  which  the  maximum  lift  coefficient  occurred,  the 
maximum  lift  coefficient  itself  did  not  increase.  The  in¬ 
crease  in  stall  angle  of  attack  is  probably,  therefore,  an¬ 
other  manifestation  of  the  effect  depressing  the  lift  curve 
slope  (Ref  19)  (c.f.  above) . 

Comparing  Fig.  8,  the  static  lift  curve,  with  Figs. 

4  through  7,  the  dynamic  lift  curves  for  the  same  test 
section  velocity,  the  following  can  be  noted.  The  most 
obvious  difference  between  the  two  types  of  lift  curves  is 
the  shape  of  the  curve  at  the  stall  point.  The  static  lift 
curve  exhibits  a  gradual  buildup  of  lift  with  increasing 
angle  of  attack  and  levels  off  as  the  stall  angle  of  attack 
is  approached.  When  the  stall  condition  is  reached,  further 
increases  in  angle  of  attack  cause  a  gradual  decline  in  lift. 
In  contrast,  the  dynamic  lift  curve  shows  a  rapid  lift  in¬ 
crease  as  angle  of  attack  increases,  with  no  indication  of 
levelling  off  as  the  stall  angle  of  attack  is  approached. 

When  the  stall  condition  is  reached,  the  loss  of  lift  is 
quite  abrupt.  Physically,  the  phenomenon  responsible  for 
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this  behavior  is  the  formation  and  growth  of  a  vortex  near 
the  leading  edge  of  the  airfoil,  followed  by  its  movement 
down  the  upper  surface  of  the  airfoil  toward  the  trailing 
edge.  Stall  occurs  abruptly  as  this  vortex  strips  off  the 
trailing  edge  into  the  freestream  (Ref  20) .  A  general  ap¬ 
preciation  for  this  effect  can  be  gained  by  comparing  the 
pressure  distribution  on  the  airfoil  immediately  before  and 
after  dynamic  stall  occurs.  Fig.  9  shows  the  pressure  dis¬ 
tribution  on  the  airfoil  just  before  dynamic  stall  occurs, 
while  Fig.  10  shows  it  a  short  time  after  dynamic  stall. 
Although  these  pressure  distributions  have  been  excerpted 
from  the  case  for  34.4  feet  per  second  test  section  velo¬ 
city  and  a  motor  voltage  of  12.00  volts,  they  are  typical 
of  those  observed  in  all  cases  in  this  study. 

Clearly,  these  pressure  distributions  for  the  dynamic 
experiment  are  different  from  those  observed  for  the  static 
case  at  the  same  angle  of  attack.  The  most  notable  feature 
of  the  pressure  distribution  for  the  dynamic  experiment  is 
the  pressure  spike  located  on  the  upper  surface  of  the  air¬ 
foil,  just  behind  the  leading  edge.  This  is  seen  to  grow 
rapidly  with  increasing  angle  of  attack  until  just  before 
stall  occurs,  when  it  undergoes  catastrophic  collapse.  A 
spike  also  occurs  for  the  static  case  at  a  comparable  angle 
of  attack,  and  ceases  to  grow  when  stall  angle  of  attack  is 
reached.  Compared  to  the  dynamic  case,  however,  the  spike 

41 


i 


Figure  9 .  Pressure  Distribution  for  Dynamic  Angle-of- Attack 
Increase  Prior  to  Dynamic  Stall 
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Figure  10.  Pressure  Distribution  for  Dynamic  Angle-of-Attack 
Increase  Following  Dynamic  Stall 


for  the  static  case  shrinks  gradually  as  angle  of  attack 
increases  beyond  stall. 

In  comparing  the  static  and  dynamic  lift  curves  fur¬ 
ther,  it  becomes  evident  that,  not  only  is  the  general  char¬ 
acter  of  the  curve  changed  in  the  neighborhood  of  the  stall 
point,  but  the  point  at  which  stall  occurs  is  changed,  as 
well.  More  precisely,  the  angle  of  attack  at  which  dynamic 
stall  occurs  is  significantly  greater  than  that  at  which 
static  stall  occurs.  This  delay  in  stall  results  in  a 
corresponding  maximum  lift  coefficient  which  is  also  signi¬ 
ficantly  greater  than  that  encountered  in  the  static  case. 
The  relationship  between  maximum  unstalled  angle  of  attack 
and  nondimens ional  angular  rate  parameter  for  the  case  of 
constant-rate  angle  of  attack  increase  was  first  determined 
by  Deekens  and  Kuebler  (Ref  3) ,  and  later  verified  by  Daley 
(Ref  10) .  Table  I  summarizes  the  results  of  the  present 
study  in  this  respect,  and  Fig.  11  shows  a  plot  of  these 
results.  Delta  stall  angle  of  attack  is  defined  as  the 
difference  between  the  dynamic  stall  angle  of  attack  and 
the  static  stall  angle  of  attack  at  the  same  Reynolds  num¬ 
ber. 

Fig.  11  shows  a  definite  disagreement  between  the  re¬ 
sults  of  the  present  study  and  those  of  Daley.  This  is,  in 
a  large  part,  due  to  the  difference  in  the  definition  of 
stall,  ^he  present  study  defines  stall  angle  of  attack  as 
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TABLE  I 


Data  Summary 


V 

00 

(Ft/Sec) 

“stall  st 
(Degrees) 

stall  dyn 
(Decrees ) 

“nd 

26.7 

16 

22.6 

.011 

26.7 

16 

25.6 

.018 

26.7 

16 

27.5 

.025 

26.7 

16 

30.7 

.032 

30.1 

16 

23.5 

.010 

30.1 

16 

24.6 

.016 

30.1 

16 

26.4 

.022 

30.1 

16 

28.7 

.029 

34.4 

16 

21.2 

.009 

34.4 

16 

24.2 

.014 

34.4 

16 

25.4 

.019 

34.4 

16 

27.7 

.025 

39.9 

17 

22.8 

.008 

39.9 

17 

25.0 

.013 

39.9 

17 

26.0 

.017 

39.9 

17 

28.6 

.022 

47.8 

17 

22.6 

.006 

47.8 

17 

23.2 

.010 

47.8 

17 

25.5 

.014 

47.8 

17 

27.1 

.018 

) 

_ 
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DELTA  STALL  ANGLE  OF  ATTACK  (DEGREES) 


0.0  0.005  0.010  0.015  0.020  0.025  0.030 
NONDIMENSIONAL  ANGULAR  RATE  «ND 


0.035 


Figure  11.  Data  Summary  Including  Daley's  Data 
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the  angle  of  attack  at  which  any  further  increase  in  angle 
of  attack  yields  no  further  increase  in  lift  coefficient. 
Daley,  however,  was  unable  to  determine  lift  coefficient, 
and  therefore  defined  flow  separation  at  the  airfoil 
quarter-chord  as  the  stall  criterion.  It  is  perhaps  signi¬ 
ficant  that  the  two  sets  of  results  exhibit  a  constantly 
increasing  divergence  with  increasing  nondimensional  angu¬ 
lar  rate.  If  Daley's  results  show  the  angle  of  attack  at 
which  quarter-chord  separation  occurs,  and  the  results  of 
the  present  study  indicate  the  actual  dynamic  stall  angle 
of  attack,  then  this  divergence  may  constitute  evidence 
that  the  length  of  time  between  quarter-chord  separation 
and  actual  stall  is  a  function  of  nondimensional  angular 
rate.  Lawrence  (Ref  15)  has  recently  completed  theoreti¬ 
cal  work  aimed  at  predicting  the  dynamic  stall  point  for 
the  case  of  constant-rate  angle  of  attack  increase.  These 
predictions  are  generally  verified  by  the  results  of  this 
experiment . 

The  final  distinction  between  the  static  and  dynamic 
stall  lift  curves  is  the  difference  in  slope.  A  clearly 
discernible  lift  curve  slope  is  not  present  in  all  twenty 
of  the  dynamic  cases,  but  in  the  cases  where  a  dynamic  lift 
curve  slope  is  recognizable,  it  seems  to  be  depressed  in 
comparison  to  the  static  case.  It  should  be  emphasized, 
however,  that  the  difference  between  the  dynamic  and  static 
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lift  curve  slopes  is  approaching  the  noise  level  of  the  ex¬ 
periment.  As  a  result,  conclusions  requiring  extensive 
quantitative  comparisons  should  probably  not  be  drawn  using 
these  results.  This  is  not  to  imply,  however,  that  general 
qualitative  trends  cannot  be  seen  here.  Tupper  (Ref  16) 
has  recently  completed  theoretical  work  intended  to  pre¬ 
dict  the  dynamic  lift  curve  slope  for  the  case  of  constant- 
rate  angle  of  attack  increase.  The  general  trend  in  the 
results  mentioned  above  seems  to  confirm  these  theoretical 
predictions . 

An  additional  investigation  was  undertaken  to  experi¬ 
mentally  determine  the  effect  of  initiating  the  airfoil  ro¬ 
tation  at  angles  of  attack  other  than  zero,  as  was  done  for 
the  main  investigation.  Airfoil  rotation  was  begun  at  -5, 

0,  5,  and  10  degrees  angle  of  attack  while  test  section  velo¬ 
city  and  motor  voltage  were  held  constant  at  29.4  feet  per 
second  and  12.00  volts,  respectively.  Since  test  section 
velocity  and  motor  voltage  were  held  constant,  the  associ¬ 
ated  nondimensional  rotation  parameter  remained  constant, 
as  well.  The  dynamic  lift  curves  obtained  in  this  investi¬ 
gation  are  presented  in  Figs.  12  through  15. 

The  variation  of  dynamic  stall  angle  of  attack  and 
maximum  lift  coefficient  seems  negligible  with  respect  to 
the  angle  of  attack  at  which  airfoil  rotation  was  begun. 

This  result  is  in  good  agreement  with  the  results  of  Law- 
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Figure  12.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =94. 6  Deg/Sec  and  ^-29  .4  Ft/Sec 

(  ®  ND=®  ’  ) 

Airfoil  Rotation  Initiated  at  a=- 5  Degrees 
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o  5  10  15  20  25 

ANGLE  OF  ATTACK  (DEGREES) 

Figure  13.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  -96.3  Deg/Sec  and  ^=29 .4  Ft/Sec 
(  «Nd=0-029  ) 

Airfoil  Rotation  Initiated  at  a=0  Degrees 
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FORCE  COEFFICIENT 
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ANGLE  OF  ATTACK  (DEGREES) 


Figure  14. 


Force  Coefficient  vs.  Angle  of  Attack 
for  0=96.1  Deg/Sec  and  ^=29 -4  Ft/Sec 

(  o  nd=®  *  ®29  ) 

Airfoil  Rotation  Initiated  at  a =5  Degrees 
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ANGLE  OF  ATTACK  (DEGREES) 


Figure  15*  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =97.1  Deg/Sec  and  Vm=29A  Ft/Sec 
(  <*nd*0.029  ) 

Airfoil  Rotation  Initiated  at  a=io  Degrees 
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rence,  which  predict  that  the  stall  point  is  a  function 
only  of  nondimensional  rotation  parameter.  No  clear  trend 
is  readily  perceptible  with  respect  to  dynamic  lift  curve 
slope ,  however . 
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VI.  Conclusions  and  Recommendations 


Conclusions 

The  present  investigation  of  dynamic  stall  has  met  all 
of  the  initial  objectives.  The  following  conclusions  can  be 
drawn  from  the  results  of  this  study. 

The  microcomputer-based  data  acquisition  system  consti¬ 
tutes  an  accurate  and  reliable  system  for  gathering  time, 
angle  of  attack,  and  pressure  information  associated  with 
the  physical  phenomenon  of  dynamic  stall.  In  addition,  the 
software  package  developed  to  automate  the  reduction  of  this 
information  is  both  efficient  and  dependable  in  accomplishing 
its  purpose. 

Based  on  the  data  obtained  for  the  NACA  0015  airfoil, 
the  following  conclusions  can  be  made  concerning  dynamic 
stall  for  nondimens ional  angular  rate  parameter  between  .006 
and  .032.  First,  there  is  a  definite  relationship  between 
the  delay  of  stall  and  the  nondimensional  angular  rate  pa¬ 
rameter,  «ND-  For  aND  =  .032,  stall  was  delayed  by  as  much 
as  14  degrees  beyond  the  static  stall  angle  of  attack  at  the 
same  Reynolds  number.  The  corresponding  aerodynamic  loads 
for  this  case  reached  C  =  -5.1  and  C.  =  1.94.  Second, 
neither  the  delay  of  stall  nor  the  corresponding  aerodynamic 
loads  appear  to  be  influenced  by  the  angle  at  which  airfoil 
rotation  is  initiated. 
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Recommendations 


The  results  of  the  present  study  show  that  there  is 
considerable  potential  for  future  dynamic  stall  work  using 
the  system  developed  for  this  experiment.  There  are,  how¬ 
ever,  some  possible  improvements,  the  implementation  of 
which  would  add  to  the  value  of  future  results. 

First,  the  stack  and  rake  now  used  in  the  smoke  tunnel 
for  smoke  trail  injection  could  be  removed.  This  might  re¬ 
duce  any  flow  irregularities  now  present  in  the  test  section, 
possibly  improving  the  quality  of  the  flow  therein. 

Second,  an  airfoil  drive  motor  capable  of  higher  angular 
rate  and  proportionately  greater  torque  output  should  be  in¬ 
stalled.  This  would  broaden  the  scope  of  the  experiment  by 
expanding  the  range  of  the  nondimensional  angular  rate. 

Third,  supplying  a  higher  excitation  voltage  to  the 
transducers  would  improve  the  signal-to-noise  ratio  as  well 
as  increase  the  percent  full-scale  deflection  at  the  output 
of  the  AIM-12  board.  This  modification  would  necessitate 
the  recalibration  of  all  transducers  at  the  new  excitation 
voltage. 

Fourth,  the  present  means  of  calibration  is  accurate, 
but  is  extremely  time-consuming.  By  constructing  an  air¬ 
tight  case  to  enclose  the  airfoil,  the  calibration  pressure 
could  be  applied  to  all  16  transducers  simultaneously. 
Coupling  with  the  proper  software  package  would  allow  auto¬ 
mation  of  the  calibration  process,  resulting  in  considerable 
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time  being  saved  for  each  calibration  run.  This,  in  turn, 
would  enable  more  frequent  calibration  of  the  transducers, 
thereby  increasing  experimental  accuracy. 

Finally,  expanded  RAM  space  in  the  computer  would  allow 
a  more  sophisticated  software  package  to  be  implemented. 
Augmented  RAM  would  also  considerably  decrease  the  time  re¬ 
quired  for  data  reduction  by  reducing  the  time  presently 
required  for  floppy  disk  input/output . 
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Appendix  A 

Tunnel  Interference  Effects 


Even  if  the  usual  Reynolds  number  comparison  criter¬ 
ion  is  met,  the  flow  around  an  airfoil  in  a  wind  tunnel  will 
generally  differ  from  that  which  would  be  observed  in  free 
air,  where  the  fluid  is  of  unlimited  extent.  The  factors 
which  produce  these  discrepancies  are  collectively  referred 
to  as  tunnel  interference  effects.  Corrections  for  these 
interference  effects  can  be  applied  to  the  parameters  meas¬ 
ured  in  the  tunnel  to  determine  the  corresponding  values 
that  would  be  obtained  for  measurement  of  the  same  para¬ 
meters  in  free  air. 

In  view  of  the  discussion  presented  by  Pankhurst  and 
Holder  (Ref  17:327-388),  three  interference  effects  were 
considered  significant  for  the  case  of  an  uncambered  air¬ 
foil  spanning  the  test  section  in  incompressible,  two- 
dimensional  flow  enclosed  by  a  test  section  having  rigid 
walls.  These  three  effects  are  two-dimensional  solid  block¬ 
age,  two-dimensional  wake  blockage,  and  two-dimensional 
streamline  curvature. 

Two-dimensional  solid  blockage  causes  an  increase  in 
axial  velocity  past  the  airfoil,  owing  to  its  partially 
blocking  the  flow  in  the  presence  of  the  wall  constraint. 
This  velocity  increment  is  a  function  of  the  chordlength, 
thickness,  and  thickness  distribution  of  the  airfoil  and  is 
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independent  of  camber  or  angle  of  attack  while  flow  remains 
attached  over  the  airfoil.  The  velocity  increment  due  to 
solid  blockage,  however,  is  much  less  than  that  which  would 
be  caused  by  direct  area  reduction.  The  reason  for  this  is 
that  the  streamlines  most  affected  are  those  at  the  tunnel 
boundary,  which  are  the  greatest  distance  away  from  the 
airfoil . 

Two-dimensional  wake  blockage  also  causes  an  in¬ 
creased  axial  velocity  past  the  airfoil,  but  does  so  for  a 
different  reason.  The  wake  behind  the  airfoil  has  a  lower 
mean  velocity  than  the  freestream  because  of  viscous  inter¬ 
action  with  the  body.  Since  the  wake  has  a  lower  mean  vel¬ 
ocity  than  the  freestream,  the  average  velocity  of  the  free¬ 
stream  must  increase  in  order  to  satisfy  continuity  in  the 
test  section. 

Two-dimensional  streamline  curvature  corrections  must 
be  applied  because  the  presence  of  the  test  section  floor 
and  ceiling  prevent  the  streamline  curvature  that  would 
normally  occur  in  the  vicinity  of  an  airfoil  in  free  air. 
Relative  to  the  flow  straightened  by  the  test  section  floor 
and  ceiling,  the  airfoil  appears  to  be  cambered,  and  there¬ 
fore  delivers  more  lift  than  it  would  at  the  same  angle  of 
attack  in  free  air. 

These  three  blockage  effects  are  quantified  mathe¬ 
matically  in  Pankhurst  and  Holder  as  follows.  The  solid 
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corrected  test  section  velocity  and  angle  of  attack  can  be 
found: 


V»CORR  V«(1  +  eSB  +  eWB) 


“  CORR  a  +  AaSC 


(13) 

(14) 


These  corrections  assume  steady-state  test  conditions 
and  were  arrived  at  using  conservative  quantities.  Tunnel 
wall  interference  effects  for  the  transient  case  are  virtually 
impossible  to  compute  accurately,  but  can  logically  be  as¬ 
sumed  smaller  than  those  for  the  static  case.  Therefore,  the 
computed  correction  factors  for  velocity  and  angle  of  attack 
define  conservative  upper  limits  on  the  error  which  would  be 
introduced  by  omitting  these  corrections .  Due  to  the  virtual 
impossibility  of  computing  tunnel  wall  interference  correc¬ 
tions  for  the  transient  case,  and  the  insignificance  of  the 
error  introduced  by  not  doing  so,  corrections  for  tunnel  wall 
interference  effects  were  not  applied  to  the  dynamic  stall 
data  gathered  in  this  study. 


Appendix  B 

Test  Section  Flow  Characterization 

The  tunnel  used  in  this  study  was  originally  intended 
for  smoke-trace  flow  visualization,  and  the  test  section  is 
only  2.75  inches  wide.  As  such,  there  was  initially  some 
concern  that  the  boundary  layers  on  the  side  walls  were  so 
thick  as  to  render  the  pressure  data  of  little  value.  A 
preliminary  investigation  was  therefore  undertaken  to  deter¬ 
mine  the  velocity  profile  of  the  transverse  boundary  layers  in 
the  test  section. 

Apparatus 

Instrumentation  used  to  measure  the  velocities  consisted 
of  the  same  Pitot-static  probe  and  Meriam  A-937  water  micro¬ 
manometer  normally  used  to  measure  test  section  velocity  in 
the  smoke  tunnel.  The  probe  was  connected  to  the  manometer 
with  two  eight-foot  lengths  of  Tygon  tubing  having  an  outside 
diameter  of  .125  inch. 

To  facilitate  placement  of  the  probe  at  the  desired  lo¬ 
cations  in  the  test  section,  a  bracket  was  constructed  from 
a  wooden  dowel  .75  inch  in  diameter.  The  dowel  was  cut  to  a 
length  of  2.625  inches  so  that  when  both  ends  were  capped 
with  a  layer  of  .125  inch  thick  double-sided  adhesive  foam 
rubber,  the  bracket  would  be  clamped  tightly  between  the 
walls  of  the  closed  test  section. 
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Three  holes,  each  .187  inch  in  diameter  and  parallel  to 
each  other,  were  bored  through  the  diameter  of  the  dowel  to 
accommodate  the  stem  of  the  Pitot-static  probe.  One  hole 
was  centered  at  the  midpoint  of  the  dowel,  and  the  remaining 
two  were  centered  at  locations  .5  and  1.0  inches  from  the 
midpoint.  Each  hole  was  equipped  with  a  set  screw  to  hold 
the  stem  of  the  probe  securely  in  the  bracket.  Thus,  each 
boundary  layer  profile  was  defined  by  four  points,  including 
the  no-slip  condition  at  the  wall  station. 

Procedure 

To  characterize  the  development  of  the  transverse 
boundary  layer  in  the  test  section,  two  bracket  locations 
were  consistently  used  in  this  investigation.  Location  A 
placed  the  tip  of  the  probe  10.0  inches  above  the  floor  of 
the  test  section  and  4.0  inches  downstream  of  where  the  test 
section  begins.  Location  B  placed  the  tip  of  the  probe  the 
same  distance  above  the  floor  as  location  A,  but  40.0  inches 
downstream  of  where  the  test  section  begins.  At  each  bracket 
location,  the  flow  velocity  was  measured  with  the  stem  of  the 
probe  positioned  successively  in  each  of  the  three  holes  in 
the  bracket. 

It  was,  of  course,  necessary  to  open  the  test  section 
each  time  the  bracket  location  or  probe  stem  position  was  to 
be  changed.  The  tunnel  motor  controls,  however,  were  set  at 
the  beginning  of  the  run  and  were  not  adjusted  until  the 
velocities  at  both  bracket  locations  and  all  three  stem  posi- 
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tions  were  measured.  Due  to  the  unusually  long  tubes  con¬ 
necting  the  probe  to  the  manometer,  care  was  taken  to  allow 
the  manometer  ample  time  to  reach  equilibrium  before  re¬ 
cording  the  reading.  It  was  also  important  to  insure  that 
the  head  of  the  probe  was  at  zero  angle  of  attack  in  order 
to  obtain  accurate  velocity  measurements.  This  was  accom¬ 
plished  by  positioning  the  probe  such  that  its  stem  was 
parallel  to  a  plumb  line  suspended  from  the  front  wall  of 
the  test  section. 

Results 

The  results  of  this  investigation  are  presented  in 
Figs.  16  through  19.  There  are  two  sets  of  velocity  pro¬ 
files,  one  for  each  test  section  velocity.  Note  that  the 
no-slip  condition  has  been  assumed  for  all  wall  points. 

From  Schlichting  (Ref  18:42),  the  approximate  thickness  of 
a  turbulent  boundary  layer  developing  on  a  flat  plate  in  the 
absence  of  a  pressure  gradient  is: 

V  1  -1/5 

f  *  °-37  (-7-)  (15> 

If  V  is  taken  as  the  velocity  at  the  centerline  of  the  test 
section  and  the  distance  down  the  test  section  is  used  for  1, 
the  predicted  thickness  of  the  boundary  layer  along  the  wall 
of  the  test  section  at  each  of  the  bracket  locations  is  as 
depicted  in  each  of  the  Figs.  As  shown,  the  predicted  and 
measured  boundary  layer  thicknesses  agree  quite  well. 


Figure  16 •  Boundary  Layer  Profile  at  Location  A 


O  Boundary  Layer  Data 
O—  — O  Boundary  Layer  Profile 

-  Predicted  Boundary  Layer  Thickness 

-  Position  of  Transducer  Closest  to  Wall 


Figure  17 .  Eoundary  La^er  Profile  at  Location  B 
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VELOCITY  ( FT/ SEC ) 


Figure  18 .  Boundary  Layer  Profile  at  Location  A 


O  Boundary  Layer  Data 
O—O  Boundary  Layer  Profile 

- —  Predicted  Boundary  Layer  Thickness 

-  Position  of  Transducer  Closest  to  Wall 


Figure  19 .  Boundary  Layer  Profile  at  Location  B 
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Conclusions 


Based  on  the  results  of  this  preliminary  investigation, 
the  effect  of  the  wall  boundary  layers  is  not  as  severe  as 
had  been  previously  supposed.  Even  at  the  lower  test  section 
velocity  and  at  the  location  downstream  of  the  airfoil  (con¬ 
ditions  which  would  give  the  thickest  boundary  layer) ,  the 
transducer  closest  to  the  wall  was  still  in  flow  having  95 
percent  of  the  velocity  at  the  centerline  of  the  test  sec¬ 
tion. 
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Appendix  C 

Transducer  Calibration 

Introduction 

Even  though  each  transducer  came  from  Endevco  complete 
with  its  own  factory  calibration,  all  16  of  the  transducers 
were  recalibrated  following  their  installation  in  the  air¬ 
foil.  As  Table  II  shows,  some  transducer  sensitivities 
changed  drastically  between  the  time  the  transducer  was 
factory  calibrated  and  the  time  it  was  first  calibrated 
after  mounting  in  the  airfoil.  Calibration  of  the  trans¬ 
ducers  was  subsequently  repeated  every  three  weeks  until 
the  experimental  portion  of  this  study  was  completed. 

After  the  initial  calibration  mentioned  above,  later  cali¬ 
brations  showed  no  transducer  undergoing  a  sensitivity 
change  of  greater  than  approximately  two  percent. 

Apparatus 

The  calibration  was  somewhat  complicated  by  the  fact 
that  the  transducers  were  embedded  in  the  airfoil,  and  the 
airfoil  could  not  be  removed  from  the  test  section  for  cali¬ 
bration.  Calibration  pressure  was  individually  applied  to 
each  of  the  transducers  using  a  small  rubber  cup  sealed  to 
the  airfoil  with  silicone  vacuum  grease.  However,  one  of 
the  central  problems  in  calibration  was  to  develop  a  method 
for  holding  the  rubber  cup  in  place  on  the  airfoil.  This 
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TABLE  II 


Transducer  Sensitivities 


Transducer 

Endevco 

Calibration 

mV 

psi 

Initial 

Calibration 

1 

173.5 

196.2 

2 

154.6 

168.4 

3 

175.6 

173.5 

4 

189.5 

226.4 

5 

187.8 

203.8 

6 

178.7 

200.1 

7 

195.7 

229.9 

8 

174.3 

208.8 

9 

160.2 

170.9 

10 

112.0 

113.9 

11 

119.6 

119.3 

12 

112.7 

112.3 

13 

131.0 

139.2 

14 

142.5 

167.6 

15 

195.6 

217.4 

16 

185.7 

217  2 
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problem  was  solved  using  the  device  pictured  in  Fig.  20. 

The  two  suction  cups  were  of  the  standard  variety  used  to 
attach  cargo  racks  to  the  roofs  of  cars,  and  had  a  nut  mold¬ 
ed  into  them  to  accept  the  bolts  shown.  The  tube  was  seal¬ 
ed  into  the  rubber  cup  using  silicone  rubber  adhesive  seal¬ 
ant.  With  the  airfoil  in  place  in  the  test  section,  the 
suction  cups  were  stuck  to  the  rear  wall  of  the  test  sec¬ 
tion,  thus  holding  the  rubber  pressure  cup  securely  in 
place  on  the  airfoil. 

The  transducers  were  calibrated  in  the  same  pressure 
range  in  which  they  operated,  the  region  of  low  magnitude, 
negative  pressures.  The  Meriam  A-937  water  micromanometer 
provided  a  means  of  simultaneously  generating  and  measur¬ 
ing  a  low  magnitude,  negative  pressure  for  the  calibration 
pressure  input  to  the  transducer.  This  was  accomplished 
by  applying  a  positive  pressure  to  the  left  leg  of  the 
micromanometer  such  that  the  water  column  in  the  right  leg 
rose  approximately  four  inches.  After  the  tube  leading  to 
the  pressure  cup  was  attached  to  the  right  leg  of  the  man¬ 
ometer,  the  pressure  to  the  left  leg  of  the  manometer  was 
released,  resulting  in  a  negative  pressure  of  approximately 
four  inches  of  water  being  applied  to  the  transducer.  To 
control  the  magnitude  of  the  pressure  applied  to  the  trans¬ 
ducer,  ambient  air  was  bled  into  the  system  through  a  tee 
fitting  midway  in  the  line  between  the  micromanometer  and 
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Figure  20.  Calibration  Device  in  Place  on  the  Airfoil 
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the  pressure  cup,  thereby  lowering  the  level  of  the  water 
column  in  the  micromanometer. 

Calibration  Procedure 

The  computer,  both  power  supplies,  and  all  three 
voltmeters  were  powered  up  and  allowed  to  warm  up  for  a 
minimum  of  four  hours  before  the  actual  calibration  was  be¬ 
gun.  When  calibration  was  ready  to  begin,  a  light  coating 
of  silicone  vacuum  grease  was  spread  on  the  lip  of  the 
pressure  cup  to  improve  the  seal  between  it  and  the  air¬ 
foil  surface,  and  the  pressure  cup  was  applied  to  the  trans¬ 
ducer  to  be  calibrated.  Next,  the  calibration  program 
l  TEST,  a  copy  of  which  is  included  in  Appendix  D,  was  exe¬ 

cuted.  With  TEST  running,  the  voltage  output  by  the  trans¬ 
ducer  was  continuously  updated  and  displayed  at  the  com¬ 
puter  terminal  in  digital  counts. 

Pressure  was  then  applied  to  the  transducer  as  prev¬ 
iously  described,  and  both  the  micromanometer  reading  and 
the  transducer  digital  output  were  then  recorded.  With  the 
calibration  system  leakproof ed,  any  fault  in  the  trans¬ 
ducer-airfoil  seal  would  be  shown  by  the  failure  of  the 
meniscus  to  reach  an  equilibrium  position.  Sufficient  air 
was  then  bled  into  the  calibration  system  to  lower  the 
water  column  approximately  one  inch,  and  the  previous  pro¬ 
cedure  was  then  repeated. 
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A  total  of  five  successively  lower  water  column 
heights  were  used  in  calibrating  each  transducer.  These 
heights  were  between  four  and  zero  inches,  spaced  at  ap¬ 
proximately  one-inch  intervals.  The  five  calibration 
points  thus  obtained  were  then  plotted  on  a  graph  having 
the  pressure  input  to  the  transducer  in  inches  of  water  on 
its  horizontal  axis  and  the  digital  pressure  reading  from 
the  transducer  in  counts  on  its  vertical  axis.  In  all 
cases,  this  calibration  curve  was  found  to  be  linear,  as 
advertised  by  Endevco.  The  slope  of  this  curve,  then, 
was  the  sensitivity  of  the  transducer  in  digital  counts  per 
inch  of  water.  The  sensitivity  in  millivolts  per  psi 
(pounds  per  square  inch)  was  calculated  using  the  follow¬ 
ing  conversion. 

mV  _  digital  counts  50  mV _  27,68  inches  of  water 

psi  -  inches  of  water  2048  digital  counts  psi  (16) 

These  transducer  sensitivities  in  millivolts  per  psi  were 
then  typed  into  the  appropriate  arrays  in  programs  RED  and 
DEN. 
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T  o  o  a  t  n  5  r  a  r  d  s  t  o  r  a  d  a  t  a  f  o  r  f  -j  r  fc  n  5  r  c  r  0  :  e  ; ;  ;  ;  r .  j 
i-ink;  UNO, STCLK. G£l  TIM.  ADIG, FGRLie/S, UNO/N/E 

IMPLICIT  INTEGER  (A-Z) 

REAL  AVSTAT < 16) , STATIC ( 16) , BAROM, TEMP, MANDM1 , MANGM2. TUNVEL 
REAL  MOTVOL , P90 , P© , RHC , DT I M , DPOSV . DPOSD . ROTRAT , VPD 
REAL  PORTU  <  1  0 ) ,  PORTL  ( 1 © ) , SENS ( 1 6 ) ,  CPU  ( 1 © ) .  CPL  (  1  © ) 

REAL  I DATAT ( 1 6 ) , NORMCO , PRESS , STI CKY 

REAL  CP ( 16 ) , AREAUT , AREALT , LNGTHU, LNGTHL, AREAU, AREAL , INTI  I, INTL 
I NTEGER  I DAT A ( 5050 ) , HOUR , CHECK , CHAN , DAY , MONTH , YEAR , X  X 
I  NTEGER  VALUE ,  CHEK ,  NS ,  N ,  A ,  D I ,  K ,  J ,  B ,  AA ,  L ,  C ,  KOUNT ,  S ,  T ,  1 1 ,  riD ,  EE ,  ' 
INTEGER  DIFANG, INK, RUNS, XXX, YYY, RRR, ZERANG, SNAP. SE'  FOT 
I NTEGER  CHECK , CHEK , CHAN , VALUE , KOUNT , Z , W , S . CCC 
I  NTEGER  1 1 ,  J.J ,  KK ,  WW ,  DD ,  X ,  V ,  Y ,  TT ,  ZZZ 

-  Lead  transducer  sensitivities  (mi  1 1 ivo 1 ts/ps i ) 

DATA  SENS/ 2® 6. 5, 174.2, 174. 9, 237. 5, 207. 205 ~2, 243. 1 

+221.2, 172.3, 113. 1, 118.3, 111.5, 13S.7, 177. 5, 221 . 9?225. 2/ 

-  Load  transducer  locations  on  upper  surface  (percent  chord) 

DAT  A  PORTU/ 0 .0,0. 0242 , 0 . 0434 , 0 . 0969 ,0.1 29 , 0 .  1  ® .  coo 

+0.388, 1 . 00®/  . 

-  Load  transducer  locations  on  lower  surface  (percent  chord ) 

DATA  PORTL/0. 0, 0.0161,0. 0319,0. 0484, 0. 0969,0. 194,0.323 

+0.686,1.00®/ 

-  Initialize  count  of  passes  to  zero. 

KOUNT =0 


Input  date,  time,  barometer,  and  room  temperature 
for  experimental  records. 

WRITE  (1,15) 

FORMAT  <■'  ENTER  DAY,  MONTH,  YEAR  SEPERATED  BY  COMMAS ■' 
READ  (1,20) DAY , MONTH, YEAR 
FORMAT  (13,13,13) 

WRITE  (1,25) 

FORMAT  < '  ENTER  TIME  (MILITARY:  XX XX  HOURS)',/) 

READ  (1,30) HOUR 
FORMAT  (15) 

WRITE  (1,35) 

FORMAT  ('  ENTER  BAROMETER  (INCHES  OF  MERCURY)'  /) 

READ  (1,40) BAROM 
FORMAT  (F7.2) 

WRITE  (1,45) 

FORMAT  ('  ENTER  ROOM  TEMPERATURE  (DEGREES  FAHRENHEIT) 
READ  (1,50) TEMP 
FORMAT  (F6.1) 

__  Reproduced  (rom 

77  ben  available  copy. 


Pate'S 


oooo  o  o  o  o 


E'  1 G 


FOR 


-  i£ r o  date.  time,  bar oms 1 3 r ,  and  roosi  tsmperat'jr; 

-  verification.  Offer-  option  -o  correct  f  au  1  fj  input , 

WRITE  < 1 , 55  > DAY . MONT H , YEAR 

55  FORMAT  (  •'  DAY:  '  ,  I  3 ,  •'  MONTH:  •'  ,13,  •'  YEAR:  •'  ,  13  ) 

WRITE  ( 1,6®) HOUR 
6®  FORMAT  (••  TIME: J, 15) 

WRITE  <1,65) BAROM 

65  FORMAT  <  •'  BAROMETER:  %  F7. 2, '  INCHES  OF  MERCURY1  > 

WRITE  (1,70) TEMP 

70  FORMAT  <•*  ROOM  TEMPERATURE:  ',  F6.  1 ,  DEGREES  FAHRENHEIT') 

WRITE  <1,75) 

75  FORMAT  <///,■'  ARE  THE  INPUTS,  ECHOED  ABOVE,  '  ■ 

WRITE  <1,30) 

80  FORMAT  <-‘  CORRECT?  IF  SO,  ENTER  A  l',/) 

READ  < 1,35) CHECK 
35  FORMAT  < 1 1  ) 

IF  < CHECK. NE. 1 )  00  TO  10 

-  Following  part  of  program  calculates  an  average  zero-input 

-  reading  for  each  transducer.  Average  is  obtained  from  10 

-  readings  of  each  transducer. 

WRITE  <1,90) 

90  FORMAT  <///,-■  THIS  PART  OF  THE  PROGRAM  OBTAINS  AVERAGE') 

WRITE  <1,95) 

95  FORMAT  (’  TRANSDUCER  ZERO- INPUT  READINGS.  WHEN  TEST-') 

WRITE  <1,100) 

100  FORMAT  <'  SECTION  VELOCITY  IS  ZERO,  HIT  RETURN  KEY') 

WRITE  <1,102) 

102  FORMAT  <■'  IN  RESPONSE  TO  "PAUSE"',///) 

PAUSE 

Initialize  all  array  elements  to  zero. 

110  CONTINUE* 

DO  120  Z= 1 , 16 
A VST AT  <  Z ) =0. 0 
120  CONTINUE 

-  Take  10  readings  from  each  transducer,  average  them  as  shown 

-  below,  then  write  these  averages  to  terminal.  Also  offer  the 

-  option  to  retake  the  average  zero-input  readings. 

DO  150  S=1 , 10 
DO  160  T=1 , 16 
CHAN=T-1 

CALL  AD  <  VALUE , CHAN , 30 ) 

AVSTAT  <  T ) =AVSTAT  <  T )  +  <  VALUE/ 10.0) 

160  CONTINUE 

150  CONTINUE 

C 
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WRITE  (1,1 55) 

155  FORMAT  <  AVERAGE  ZERO- INPUT  READINGS  ROLLON',/' 

DO  130  W=  1 ,  16 

WRITE  <1, 165)W, AVSTAT(W) 

165  FORMAT  ( •'  TRANSDUCER '  ,  13,  J  AVERAGE  STATIC  READING:  '  ,  F6.  0 ) 

180  CONTINUE 

WRITE  (1,177) 

177  FORMAT  <///,-'  TO  PROCEED  WITH  THE  PROGRAM,  ENTER  A  1',/) 

READ  (1,178) XX 
173  FORMAT  (12) 

IF  (XX.NE. 1 )  GO  TO  110 


C 

c 


Enter  manometer  reading,  motor  voltage,  and  9®  and  ® 
degree  angle  of  attack  voltages  for  experimental  records. 
Test-section  velocity  is  also  computed  as  shown  below. 


WRITE  (1,185) 

135  FORMAT  (//////,' *************************** *NOW  TURN  ON  THE 

+  TUNNEL****-*"*1-************  **•»•#•*•*•*■»■*  •',//////) 

137  WRITE  (1,190) 

190  FORMAT  (-•  ENTER  ROOM  PRESS.  MINUS  TUNNEL  ST  AT.  PRESS. 

+  (INCHES  OF  WATER)',/) 

READ  (1, 1931MANOM1 

195  FORMAT  (F3.4) 

WRITE  (1,200) 

200  FORMAT  (•'  ENTER  TUNNEL  TOTAL  PRESS.  MINUS  TUNNEL  STATIC  PRE 
+  (INCHES  OF  WATER)',/) 

READ  ( 1 , 195)MAN0M2 

205  FORMAT  (F3.4) 

WRITE  (1,210) 

210  FORMAT  (•'  ENTER  MOTOR  VOLTAGE  (VOLTS)',/) 

READ  ( 1 ,215) MOT VOL 

215  FORMAT  (F6.2) 

WRITE  (1,220) 

220  FORMAT  ( •’  ENTER  90  AND  0  DEGREE  VOLTAGES,  RESPECTIVELY'’,/) 
READ  ( 1,225 )P90,P0 
225  FORMAT  (2F7.4) 

RHO- ( BAROM*70 . 45  > / ( 1 7 1 6 . 0* ( 460 . 0+TEMP ) ) 

TUNVEL“SQRT ( ( 2 . 0* ( 5 . 204*MANOM2 ) ) /RHO ) 


Echo  manometer  readings,  tunnel  velocity,  motor  voltage  and 
90  and  0  degree  angle  of  attack  voltages  for  verification, 
offer  option  to  correct  faulty  input. 

WRITE  ( 1,230 1MANQM1 
230  FORMAT  (  '  MANOMETER  ONE:  ',F8.4,-'  INCHES  OF  WATER') 

WRITE  ( 1 , 233 ) MAN0M2 

233  FORMAT  (  '  MANOMETER  TWO:  ',F3.4,'  INCHES  OF  WATER'') 
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WRITE  ( 1 , 235; TUNVEL 

223  -ORMAT  •  TUNNEL  VtLOCITY:  '.F7.2,'  ,-T/SEC'  ) 

WRITE  ( 1,2^0 )MGTVOL 

240  FORMAT  ('  MOTOR  VOLTAGE:  1-6.2,-'  VOLTS') 

WRIIE  (1,245)F'90,F0 

245  FORMAT  (•'  P  90:  •'  ,  F7.4,-'  VOLTS  P0:  ',F7.4,'  VOLTS') 

WRITE  <1,75) 

WRITE  (1,30) 

READ  < 1 ,85)CHEK 
IF  (CHEK.NE.l)  GO  TO  137 

-  Initialize  number  or  runs  to  zero,  and  then  increment  this 

-  number  by  one  each  run  thereafter. 

RUNS=0 

>50  CONTINUE 

RUNS=RUNS+1 
>55  CONTINUE 


WRITE  <1,257) RUNS 

>57  FORMAT  (////,-'  *-*******RETURN  AIRFOIL  TO  ZERO  ANGLE  Qe 
+  ATTACK  IN  PREPARATION  FOR  RUN-'  ,  12,  '  *********  , //// ) 

NS=0 
KOUNT =0 
WRITE  (1,260) 

260  FORMAT  ( '  ENTER  NUMBER  OF  SAMPLES  (MULTIPLE  OK  13, 

+  5040  MAXIMUM)  •'  ,  /  ) 

READ  (1,265) NS 
265  FORMAT  (15) 

WRITE  (1,270) NS 

270  FORMAT  <//,'  '  , 25X, ' NS: ' , 15, // ) 

-  In  the  next  segement,  the  operator  is  given  the  choice 

-  between  manual  and  automatic  trigger. 

WRITE  (1,273) 

273  FORMAT  ( J  DO  YOU  WANT  MANUAL  OR  AUTOMATIC  TRIGGER? 

+  ( l=AUTO,  2=MANUAL )' ,/) 

READ  (1,277) SELECT 
277  FORMAT  (12) 

IF  < SELECT. NE. 1 )  GO  TO  282 


The  program  segement  below  is  the  automatic  trigger. 
The  program  stays  in  the  280  loop  below  until  ZERANG 
and  VALUE  differ  by  2  or  more  digital  counts. 

When  this  occurs,  due  to  rotation  of  the  airfoil,  the 
program  continues  on  to  line  number  235. 

CALL  AD (VALUE, 0,84) 

ZERANG=VALUE 
230  CALL  AD (VALUE, 0, 34) 

SNAP= I ABS  <  VALUE-ZERANG ) 
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B I  G 


FOR 


IF  ■:  SNAP.  LE.  1 )  00  ro  280 
IF  <  SELECT . EG). 1 )  GO  70  03  5 
282  CONTINUE 
PAUSE 

-  STCLK,  below,  will  count  up  to  32,768  time  clicks,  each  click 

-  being  .0010046  seconds  long.  Therefore,  STCLK  can  only  time 

-  an  event  that  lasts  for  no  more  than  about  32  seconds. 

285  CALL  STCLK 

-  The  following  part  of  the  program  reads  and  stores  the  time 

-  obtained  from  subroutine  GETTIM,  as  well  as  oosition  and 

-  pressure  information  obtained  from  the  potentiometer  and 

-  pressure  transducers,  respectively.  This  position  and  precsu 

-  information  is  obtained  through  subroutine  ADIO. 

WRITE! 1 , 290) 

290  FORMAT  (///,'  '  ,  20  X ,  •*  START  I NG  TO  TAT  E  DATA  ',///) 

DO  320  J= 1 , NS , 1 8 
KGUNT =KOUNT.+ 1 
CALL  GETTIM (TIME) 

IDATA( J)=TIME 
CHAN=0 

CALL  AD (VALUE, CHAN, S4) 

I DATA ( J+l) -VALUE 
DO  300  K=1 , 16 
CHAN=K- 1 

CALL  AD ( VALUE , CHAN , S0 ) 

DI=K+J+1 
IDATA(DI )=VALUE 

-  The  do  305  loop  below  is  a  simple  delay  loop  and  can  be  used 

-  govern  the  number  of  time  clicks  that  will  elapse  between  pas 

-  One  click  =  .0010046  seconds  and  one  pass  is  defined  as  one 

-  complete  cycle  around  the  18  data  channels.  To  find  what  ZZ 

-  should  approximately  be,  use  the  following  equation: 

CLICKS/CYCLE=(0.46)*ZZ  +  A. 0 

-  CAUTION:  Making  ZZ  too  large  will  result  in  data  points  that 

-  are  too  far  apart,  thereby  invalidating  the  linear  interpo- 

-  lation  that  is  done  in  the  data  reduction  program,  RED. FOR 

DO  305  ZZ=1 , 40 
C305  CONTINUE 

300  CONTINUE 

320  CONTINUE 

WRITE  (1,330) RUNS 

330  FORMAT  ( •*  ‘,15X,'DATA  GATHERING  COMPLETE  FOR  RUN',  12,//) 

WRITE  ( 1 , 340 ) KOUNT 

340  FORMAT  ( •’  NUMBER  OF  PASSES  =  ',16,//) 
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FOR 


N=KOU.MT  *  1 S 
WRITE  (  1 ,  343  )  N 

343  FORMAT  ('  NUMBER  OF  I  DATA  ELEMENTS”  '  .  16,  /  •'  > 

VPD”  <  P90-P0 )  /90.  0 

DTIM”< I DATA ( 1801 ) -IDATA ( 343 ) ) * ( O. O01OO46 > 

DF'OSV=<  <  IDATA(  1802) -IDATA (344)  > /4096. 0)*10.  0 
DPOSD=DPOS V / VPD 
ROTRAT =DPOSD / DT I M 
WRITE  ( 1,410) ROTRAT 

410  FORMAT  <  AIRFOIL  AVERAGE  ROTATION  RATE:  '  ,  F5.  2,  J  DEG/SEC'  , /// 

-  Options  are  now  offered  to  list  the  IDATA  array  at  the 

-  terminal,  to  write  this  array  to  disk,  and  to  reoeat  the 

-  data  run. 

WRITE (1,345) 

345  FORMAT  (■•  DO  YOU  WANT  TO  LIST  THE  IDATA  ARRAY? (  Y~1  )*,//) 

READ < 1 ,347)AA 
347  FORMAT  (12) 

IF  (AA.NE.l)OO  TO  350 
DO  420  XXX=130,N, ISO 
YYY=XXX-179 

WRITE  ( 1 , 360 ) ( I DATA (L ) , L=YYY , X X X ) 

360  FORMAT  <917 ) 

PAUSE 

120  CONTINUE 

WRITE  (1,351) 

351  FORMAT  <//) 

350  WRITE (1,355) 

355  FORMAT ( '  DO  YOU  WANT  TO  WRITE  TO  DISK?( Y=1 ) ' , //  ) 

READ  <1,347)8 
IF  (B.EQ.l)  GO  TO  390 
WRITE  (1,375) RUNS 

375  FORMAT  ('  DO  YOU  WANT  TO  REPEAT  RUN ’,12,'?  (Y”l  )',//) 

READ  ( 1 , 3S0 ) C 
380  FORMAT  (12) 

IF  (C.EQ.l)  GO  TO  255 

IF  (C.NE.l)  GO  TO  4800 

390  CONTINUE 

-  The  following  part  of  the  program  writes  pertinent 

-  information  to  file  RAWDATAODAT  on  disk. 

IF  (RUNS. NE. 1 )  GO  TO  705 

CALL  OPEN  <3, 'RAWDATAODAT' ,2) 

WRITE  (3,500) 

500  FORMAT  ('  DAY' , 10X , ' MONTH' , 9X, ' YEAR' , 9X , ' TIME' ) 

WRITE  (3,510) DAY , MONTH , YEAR , HOUR 
510  FORMAT  < 13,  1 1 X, 13, 1 1 X, 13, 9X, 15, / ) 

WRITE  (3,520) 

520  FORMAT  ('  TEMPERATURE' , 14X , ' BAROMETER' ) 

WRITE  (3,530) TEMP , EAROM 
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J0  S*  ORiiA  !  (  2/  ,  F  i- .  1  ,  1  ttX  ,  r  7  .  2 ,  /  i 

wR I  i  £  i'3,  540  5 

540  FORMAT  ( '  MANOMETER  1 ’ , 22X . ' MANOME I ER  2  '  > 

WRITE  ( 3 . 545 ) MANOM 1 , MAN0M2 
545  FORMAT  ( 2X , F8. 4, 25X, F8. 4, / ) 

WRITE  (3,550) 

550  FORMAT  ( '  TUNNEL  VELOCITY-1 ,22X,  'MOTOR  VOLTAGE') 

WRITE  ( 3 , 555 ) TUNVEL , MOTVOL 
555  FORMAT  ( 4X , F7. 2, 31 X , F6. 2, / ) 

WRITE  (3,560) 

560  FORMAT  ( •'  90  DEG.  VOLTAGE ' ,  16X,  '  0  DEG.  VOLTAGE') 

WRITE  (3, 570 )P90, P0 
570  FORMAT  <5X, F7. 4, 23X , F7. 4, / > 

WRITE  (3,530) 

530  FORMAT  < '  NUMBER  OF  PASSES' , 10X ,' NUMBER  OF  IDATA  ELEMENTS') 
WRITE  (3,590) 

590  FORMAT  ( 5X , ' ( KOUNT )  •’  ,  26X ,  '  <  N ) '  > 

WR I TE  ( 3 . 6@0 ) KOUNT , N 
600  FORMAT  ( 3X ,16, 26X ,16,//) 

WRITE  (3,610) 

610  FORMAT  <•'  AVERAGE  ZERO- INPUT  READINGS  GIVEN  BELOW',/) 

WRITE  (3, 620) AVSTAT ( 1 ) , AVSTAT (2 ) , AVSTAT ( 3 ) , AVSTAT ( 4 > 

WRITE  (3, 620) AVSTAT (5) , AVSTAT (6) , AVSTAT (7) , AVSTAT (S) 

WRITE  (3,620) AVSTAT ( 9 ) , AVSTAT ( 1 0 ) , AVSTAT < 1 1 ) , AVST AT ( 1 2 ) 
WRITE  < 3 , 620 > AVSTAT (13), AVSTAT (14), AVSTAT (15), AVSTAT ( 1 6 ) 

620  FORMA'  ‘ c9. 3, 5X , F9. 3, 5X , F9. 3, 5X , F9. 3 ) 

WRITE  3,660) 

660  FORMAT  (///) 

705  CONTINUE 

-  The  part  of  the  program  below  writes  the  collected  data 

-  to  disk,  in  unformatted  form,  under  the  filename 

-  RAWDATA1DAT,  RAWDATA2DAT,  .  .  .  .  ,  RAWDATA5DAT,  depending 

-  on  the  value  of  the  variable  RUNS.  To  view  the  data  files 

-  that  are  in  unformatted  form,  use  program  LOOK. 

IF  ( RUNS . EQ . 1 )  GO  TO  710 
IF  ( RUNS. EQ. 2 )  GO  TO  720 
IF  (RUNS. EG. 3)  GO  TO  730 
IF  (RUNS. EQ. 4 )  GO  TO  740 
IF  ( RUNS . EQ . 5 )  GO  TO  750 

710  CONTINUE 

CALL  OPEN  <  4 , ' RAWDAT A 1 DAT ' , 2 ) 

WRITE  ( 4 ) ( IDATA ( L ) , L= 1 , N ) 

GO  TO  760 
720  CONTINUE 

CALL  OPEN  (5, ' RAWDATA2DAT' ,2) 

WRITE  <5)<IDATA(L),L=1,N) 

GO  TO  760 
730  CONTINUE 

CALL  OPEN  ( 6 , ' RAWDAT A3DAT '  ,  2  > 
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UK  I TE  <  6  '  (  I  BATA  i  L  )  .  L~  t  ,  M  > 

GO*  TO  76U 
740  CONTINUE 

CALL  OPEN  ( 7 ,  '  RAWDA TA40AT  '  ,2) 

WRIIE  ( 7 ) < I DATA ( L ) , L~ 1 , N ) 

GO  TCi  760 
750  CONTINUE 

CALL  OPEN  (8, ' RAWDAT A5DAT 1 , 2 ) 

WRITE  ( 8 ) < I BAT A ( L ) , L= 1 , N  > 

GO  TO  760 
760  CONTINUE 

IF  < RUNS . NE . 5 >  GO  TO  250 
WRITE  (1,2345) 

2345  FORMAT  <//////,•'  FOLLOWING  PART  OF  PROGRAM  GIVES  STATIC 
+  NORMAL  COEFF.  FOR  STATIC  ALPHA',//////) 

C 

C  -  The  remaining  portion  of  the  program  takes  and  processes 

C  -  data  for  static  angle  of  attack  lift-curves. 

C 

240©  CONTINUE 

WRITE  (1,2450) 

2450  FORMAT  <  •’  ENTER  NS  (MULTIPLE  OF  18,  LESS  THAN  OR 
+  EQUAL  TO  396  )•',/) 

READ  < 1 , 2150 ) NS 
2150  FORMAT  (14) 

KOUNT =0 
WRITE  (1,2000) 

2000  FORMAT  (////,•'  HIT  RETURN  TO  START  DATA  COLLECTION',/) 
PAUSE 


STCLK,  he  1  w,  will  count  up  to  32,768  time  clicks,  each  click 
being  .0010046  seconds  long.  Therefore,  STCLK  can  only  time 
an  event  that  lasts  for  no  more  than  about  32  seconds. 

CALL  STCLK 

WRITE ( 1,2100) 

2 1 00  FORMAT  (///,-'  '  ,  20X ,  '  START  I NG  TO  TAKE  DATA  ',///) 

DO  2200  J=1 , NS, 18 
KOUNT =KOUNT  + 1 
CALL  GETTIM(TIME) 

IDATA(J)=TIME 

CHAN=0 

CALL  AD (VALUE, CHAN, 84) 

I DAT A ( J+ 1 ) =VALUE 
DO  2300  K=1 , 16 
CHAN=K- 1 

CALL  AD ( VALUE, CHAN, 80 ) 

DI=K+J+1 
I DAT A ( D I ) =V ALUE 
2300  CONTINUE 
2200  CONTINUE 
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N=K0UNT*18 
UR  lit  (1,2  j  -  >■:> >  N 

2500  FORMAT  ('  NUMBER  OF  I  DATA  ELEMENTS”  ‘,16,/.'' 

PAUSE 

C 

C  -  Time-average  data 

C 

DO  2550  S=l,16 
IDATAT ( S ) =0. 0 
2550  CONTINUE 

DO  2600  1 1  =  1, N, 18 
DO  2700  JJ=3,18 
TT=II+.JJ 

I  DAT  AT  ( J  J-2 )  =  (  <  IDATA(TT-1  )  ) /(COUNT  )  + IDATAT  ( J.J-2 ) 

2700  CONTINUE 
2600  CONTINUE 
C 

C -  Compute  the  pressure  c o efficients 

C 

DO  2300  KK= 1 , 16 

ST I CKY=AVSTAT ( KK ) -2043 . 0 

PRESS*  (  (  ( I  DAT  AT  <  KK )  -ST  I CKY )  -2048 . 0 )  / 2048 .  0 )  *  ( 50 . 0 /  Sir  NS:  ( KK )  ) 

CP ( KK ) ~ ( PRESS+ ( MANOM 1 /27 . 68 ) ) / ( MAN0M2/27 . 63 ) 

2300  CONTINUE 
C 

C  -  The  next  loop  defines  the  pressure  distribution  on  the  upper 

C  -  surface  of  the  airfoil,  leading  edge  to  trailing  edge. 

C  -  Pressure  coefficient  is  assumed  to  be  zero  at  the  trailing  edge 

C 

WRITE  (1,2900) 

2900  FORMAT  ('  UPPER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L.E.  TO  T.E.,  ARE  GIVEN  BELOW’,) 

DO  3000  V=1 , 9 
CPU(V)=CP(V) 

3000  CONTINUE 

CPU< 10) =0.@ 

DO  3100  V=l,10 
WRITE  ( 1 , 3200 ) V , CPU < V ) 

3200  FORMAT  ( '  CPU ' ,  13,  •'  =  ■*  ,  F8.4) 

3100  CONTINUE 
C 

WRITE  (1,3300) 

3300  FORMAT  (/,*  LOWER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L.E.  TO  T.E.,  ARE  GIVEN  BELOW' ) 

CPL ( 1 )=CP( 1 ) 

DO  3400  W=2, 8 

DD=13-W 

CPL< W)=CP(DD) 

3400  CONTINUE 

CPL ( 9 ) =CPU ( 10) 

DO  3500  W=1 , 9 

WRITE  ( 1 , 3600 ) W , CPL ( W ) 
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Oil  O 
3500 


3300 


3900 

H 

4000 

3700 


FukMA  1  <  ■' 
CO  NT  1  MUE 


CF’L '  , 


i  ri  e  folio  to  i  n  g  1  o  op  l  -i  t  e  g  r  a  t  a  s  t  h  e  u  p  p  a  r  p  r  a  c  <j  r  a 

distribution  using  the  trapezoidal  rule. 

AREAUT=0.0 
DO  3700  X=l,9 
LNGTHU=PGRTU (  X  + 1  > -PORTU (  X ) 

IF  ( (ABS(CPU(X+1)-CPU(X) > >.GT. (ABS( (0.01 >*CPU(X> ) > )  GO  TO  3800 
AREAU= ( 0 . 5 ) * ( CPU ( X  + 1 ) +CPU <  X ) ) *LNGTHU 

IF  ( (ABS(CPU(X+1)-CPU(X> > >.LE. <ABS< (0.01 )*CPU(X) > ) )  GO  TO  4000 
I NTU= ( PORTU ( X ) -PORTU <  X  + 1  >  ) *CPU ( X ) / ( CPU ( X  + 1 ) -CPU ( X  >  > 

IF  ( INTU.LT.LNGTHU)  GO  TO  3900 
AREAU*  ( .  5 )  *  ( CPU  (  X  +  l )  +CPU  ( X  )  >  *LNGTHU 
IF  (  ( INTIJ) .  GE.  <  LNGTHU )  )  GO  'TO  4000 
AREAU= ( ( . 5 ) * I NTU*CPU ( X ) ) + 

(  <  .  5 )  *  ( LNGTHU- 1  NTI.I )  *CPU  (  X+ 1  )  > 

AREAUT =AREAUT +AREAU 
CONTINUE 


-  The  following  loop  integrates  the  lower  pressure 

-  distribution  using  the  trapezoidal  rule. 

AREALT=0. 0 
DO  4100  Y=l,8 
LNGTHL=PORTL ( Y+ 1 ) -PORTL ( Y ) 

IF  ( (ABS(CPL(Y+1)-CPL(Y) ) >.GT. (ABS( (0.01 )*CPL(Y) ) ) )  GO  TO  4200 
AREAL= < . 5 ) * ( CPL ( Y+ 1 ) +CPL  < Y) ) *LNGTHL 

IF  ( ( ABS ( CPL ( Y+ 1 ) -CPL ( Y ) ) ) . LE . ( ABS ( (0.01) *CPL ( Y ) ) ) >  GO  TO  4400 
4200  I NTL= ( PORTL ( Y ) -PORTL ( Y+ 1 > ) *CPL ( Y  > / ( CPL ( Y+ 1 ) -CPL  <  Y  >  ) 

IF  ( < I NTL ) . LT . ( LNGTHL ) )  GO  TO  4300 
AREAL* ( . 5 ) * ( CPL ( Y+ 1 ) +CPL ( Y) ) *LNGTHL 
IF  ( (INTL) .GE. (LNGTHL) )  GO  TO  4400 
1300  AREAL* ( ( . 5 ) * I NTL*CPL ( Y ) ) + 

+  < ( . 5 ) * ( LNGTHL- I NTL ) *CPL ( Y+ 1 ) ) 

1400  AREALT *AREALT +AREAL 
1 1 00  CONTINUE 


4300 

H 

4400 

4100 


NORMCO=AREALT-AREAUT 
WRITE  ( 1,4500 )NORMCO 

4500  FORMAT  (/,'  NORMAL  FORCE  COEFFICIENT*' , F3. 5, / ) 

PAUSE 

-  Option  now  offered  to  repeat  the  data  run. 

WRITE  (1,4600) 

4600  FORMAT  ('  DO  YOU  WANT  TO  REPEAT  THE  RUN?  (N*l  >•’,/) 
READ  ( 1 , 4700 ) CCC 
4700  FORMAT  (12) 

IF  (CCC. NE. 1 )  GO  TO  2400 
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I NTEGER  I  DATA  (  5050  > ,  I  DATA  1  C  3« )  ,  i  DA"7 A 2  <  300  >  .  PASS  ■'  30 
I  NT  EGER  fi,R.  X,  Y ,  V ,  W ,  S .  I .  .J ,  L .  AA ,  PR .  00 

INTEGER  RR, SS, TT, UU, VV, WW, XX, YY, ZZ, RUN, TRAP, PAZZ, DIV, KIJMEL 
I NTEGER  ELEM 1 , ELEM2 , DAY , MONTH , YEAR , HOUR , CHANG 1 , CHANG2 
INTEGER  DD , EE , FF , HH , LL , NN 

REAL  PORTU  < 1 0 ) , PORTL (10), CP (16), CPU (10), CPL (10), SENS (16) 
REAL  PRESS (16), REDAT ( 40 ) , P90 , P0 , TEMP , BARQM , MANOM 1 , MAN0M2 
REAL  TUNVEL , MOTVQL , AVST AT (16), AREAUT , AREALT , RE , RHO , MU 
REAL  VPD, AOA, INTU, INTL, NORMCO, TUNG, LNGTHU, LNGTHL 
REAL  AREAU , AREAL , DT I M , DPOSD , DPOSV , ROTR AT , NDRATE 
REAL  REDATC ( 40 ) , I DAT AT ( 1 1 00 ) 

- Load  transducer  sensitivities  (mi  1 1 i vo 1 ts/psi ) 

DATA  SENS/206.5,  174.2, 174.9,237.5,207.0.205.2,243. 1. 

+221.2,  172.  S,  113.  1,  US.  3,  111.5,  13S.7,  177.5,221.9,225.2/ 

- Load  transducer  locations  on  upper  surface  (percent  chord  ‘ 

DATA  PORTU/0. 0, 0. 0242, 0. 04S4, 0.0969, 0. 12°, 0. 194. 0.323. 0. !0t 
+0. 3 S3, 1 . 000/ 

-  Load  transducer  locations  on  lower  surface  (percent  chord) 

DATA  PORTL/0. 0, 0.0161, 0.0319,0.0484,0.0969,0.  194 , 0.  323, 
+0.636,1.000/ 

WRITE  (1,5) 

5  FORMAT  (///,-•  *****THE  DATA  FILES  TO  BE  REDUCED  MUST  EE  CN 
+  DISK  DRIVE  B  AND  MUST  BE  NAMED*****"*  ) 

WRITE  (1,6) 

6  FORMAT  (’  ********** ***RAWDATA0. DAT,  RAWDATA1.DAT, . 

+  RAWDATA5 .  DAT*************  •*,///> 

-  Read  raw  data  from  RAWDATAODAT  on  drive  B. 

CALL  OPEN  ( 3 ,  RAWBATA0DAT  '  ,  2 ) 

READ  (3,10) DAY . MONTH , YEAR , HOUR 
10  FORMAT  (/,  13,  11X,  13,  UX,  13, 9X,  15) 

READ  ( 3, 20 ) TEMP, BAROM 
20  FORMAT  (//,2X,F6. 1, 18X.F7.2) 

READ  ( 3 , 30 ) MANOM 1 , MAN0M2 
30  FORMAT  (//,2X,F3.4,25X,F8.4) 

READ  (3,40) TUNVEL , MOTVOL 
40  FORMAT  (//,4X,F7.2,31X,F6.2) 

READ  <3,50>P90,P0 
50  FORMAT  ( // , 5X , F7. 4, 23X , F7. 4 ) 

READ  (3,60) KOUNT , N 
60  FORMAT  ( ///, 3X, 16, 26X, 16) 

READ  ( 3 , 70 ) AVSTAT ( 1 ) , AVSTAT ( 2 ) , AVSTAT ( 3 ) , AVST AT ( 4 ) 

READ  (3,75) AVSTAT ( 5 ) , AVSTAT  <  6 ) , AVSTAT ( 7 ) «  AVSTAT ( 8 ) 

READ  ( 3, 75 ) AVSTAT ( 9 ) , AVSTAT (10), AVSTAT (11), AVSTAT (12) 

READ  (3,75) AVSTAT (13), AVSTAT ( 1 4  > , AVSTAT (15), AVSTAT ( 1 6 > 

70  FORMAT  ( / / / / , F9 . 3 , 5X , F9 . 3 , 5X , F9 . 3 , 5X , F9 . 3 ) 

75  FORMAT  ( F9. 3, 5X , F9. 3, 5X , F9. 3, 5X , F9. 3 ) 
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C  -  The  to  *  :  owing  par  t  or  the  c  ro  gram  r  r- g :  -  s-  •  -  •  >  j 

C - attack  f  r  om  cata  sett  i,  11  21,  .  .  .  ?*■:..  r 

C - IDATA  array  (where  a  Cata  set  is  defir.ec  as  the  13  data 

C  -  elements  read  in  one  sampling  pass).  These  data  sets  were 

C  -  previously  stored  in  the  IDATA  array  in  the  first  data  file, 

C  -  RAWDATA1DAT,  by  the  data  gathering  program  BIG.  After  these 

C  -  angles  of  attack  are  read  from  the  IDATA  array  in  RAWDATADAT1, 

C  -  they  are  stored  sequentially  in  the  array  IDATA1. 

C 

ELEM1=0 

CALL  OPEN  ( 4 ,  J  RAWDAT A 1  DAT  ,  2 ) 

READ  <  4 ) ( I DATA (L ) , L~1,N) 

DO  500  PP=2,N,180 
ELEM 1 ~ELEM 1+1 
I  DATA  1  (ELE'il  )=IDA'TA(PP) 

5 0y  CONTINUE 

C 

C  -  The  DO  65@  loop  below  begins  by  branching  out  to  one  of  4 

C  -  program  locations  depending  on  the  number  of  times  the 

0  -  computer  has  been  through  the  650  loop.  The  program  segemsnts 

C  -  beginning  at  510,525,535  and  545  all  perform  the  same 

C  -  operation,  but  do  so  on  different  data  files,  in  turn. 

-  The  operation  performed  in  all  4  segements  is  as  follows. 

-  The  loop  contained  in  each  segement  reads  tne  raw  digital 

-  angle  of  attack  from  every  data  set  in  the  array  IDATA. 

-  These  data  sets  were  previously  stored  in  the  second  through 

-  fifth  data  files,  RAWDATA2DAT  through  RAWDATA5DAT  by  the  data 

-  gathering  program  UNO.  After  these  angles  of  attack  are  read 

-  from  the  IDATA  array  in  the  proper  data  file,  they  are  stored 

-  in  the  array  IDATA2.  Note  that  the  contents  of  the  1DATA2  array 

-  are  replaced  each  time  the  DO  650  loop  is  travelled. 

DO  650  RUN*2,5 
ELEM2=0 

IF  (RUN. EQ.2)  GO  TO  510 
IF  (RUN. EG. 3)  GO  TO  525 
IF  (RUN. EG. 4)  GO  TO  535 
IF  (RUN. EQ. 5)  GO  TO  545 
510  CONTINUE 

CALL  OPEN (5, 1 RAWDATA2DAT %  2  > 

READ ( 5 ) ( I DAT A  <  L ) , L= 1 , N ) 

DO  520  QQ=2, N, 13 

ELEM2=ELEM2+1 

I DATA2 ( ELEM2 ) = I DATA ( GO ) 

520  CONTINUE 
GO  TO  550 
525  CONTINUE 

CALL  OPEN  ( 6 ,  RAWDATA3DAT '  ,  2 ) 

READ ( 6 ) ( IDATA(L),L*1,N) 

DO  530  QQ=2,N, 18 
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ELEM2-ELEH2+  i 
i  D .  Ai  ■  ELElv12  )  ~  I  DA  i  A  (  yij  ) 

5  30  COM  r I NUE 

GO  TO  350 

535  CONTINUE 

CALL  OPEN (7, * RAWDAT A4DAT ‘,2) 

READ ( 7 ) ( IDATA(L) ,L=1,N) 

DO  540  Q0=2, N, 18 
ELEM2=ELEM2+ 1 
I DATA2 ( ELEM2 ) = I DATA ( QQ ) 

540  CONTINUE 
GO  TO  550 
545  CONTINUE 

CALL  OPEN <8.  ' RAWDAT A5E AT ■*  ,  2 ) 

READ ( S ) ( IDATA<  L ) , L=1 , N  > 

DO  550  QQ-2, N,  13 
ELEN2~ELEM2+ 1 
IDATA2 < ELEM2 ) = I DATA ( QO > 

550  CONTINUE 
C 

C  -  We  now  have  an  array,  IDATA1,  containing  every  tentn  angle 

C  -  attack  element  from  disk  file  RAWCATA1DA T.  We  also  have  an 

C  -  array,  IDATA2,  containing  every  angle  of  attack  element  frcn 

C  -  one  of  the  four  remaining  disk  files.  The  program  segement 

C  -  below  compares  each  element  in  10ATA1  with  every  element  in 

C  -  IDATA2,  in  turn.  When  the  closest  match  is  found,  this  info- 

C - mat  ion  is  stored  in  the  PASS  array  for  use  in  the  ne;ct  oro- 

C - qram  seqement. 

C 

RR=0 

DO  560  SS=1 , ELEM1 

CHANG2=4096 

TRAP=0 

DO  570  TT=1 , ELEM2 
RR=RR+1 

CHANG 1 * I AB3 ( I DATA 1 ( SS ) - I DATA2 ( TT ) ) 

IF  < CHANG 1.GE.CHANG2)  GO  TO  580 
CHANG2=CHANG 1 
GO  TO  570 
580  TRAP=TRAP+1 

IF  (TRAP. NE. 1 )  GO  TO  570 
PASS(SS)«RR-1 
570  CONTINUE 
560  CONTINUE 

-  The  next  part,  up  to  and  including  line  610,  reads 

-  two  data  sets  every  ten  data  sets,  ie,  1,2,  11,12.  21,22, 

-  etc.,  from  disk  file  RAWDATA1DAT,  and  then  repacks  each 

-  data  set  sequentially  in  array  IDATAT,  ie,  data  sets 

-  1,2,11,12,21,22,  ....  in  IDATA  become  data  sets 

-  1,2, 3, 4, 5, 6,  .  .  .  .  in  IDATAT. 
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Ktt-iD  1  'i  )  (  I  DA  i  h  \  L  )  ,  L  - 1 ,  N ) 

CAi_L  OPEN  (  9 ,  ’  AVRuDAT ADAT  '  ,  A  ; 
PAZZ=® 

DO  61®  UU= 1 ,  N ,  1 80 

PAZZ=PAZZ+1 

DO  62®  W-l ,  36 

DD=(UU+W)-1 

WW= W+  <  36*  ( PAZZ- 1 )  ) 

1 DATAT  <  WW  >  =  I DATA ( DD ) 

620  CONTINUE 
610  CONTINUE 

NN“PAZZ*36 

WRITE ( 9 ) ( IDATAT(LL) ,LL=1,NN) 
630  CONTINUE 


C 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 


The  foil  o w i  ri g  p r  o g r a m  3 e  g  e m en!:;  beginning  o n 
69®,  and  720  perform  the  same  operation,  bur 
different  disk  files,  in  turn.  The  operation 
follows.  The  information  previously  obtained 


1 l nes  6-0 . 665 , 

O  '0  s  o  0  fi 
performed  is  as 
and  s  t  o r e d  i  n 


the  array  PASS  is  now  used  to  determine  which  data  sets  will 
t>e  read  (from  whichever  disk  file  of  the  remaining  four  that 
is  then  being  worked  on),  then  added,  element  by  element,  to 
the  element  totals  in  the  array  I DATAT.  Note  that  the  contents 
of  the  array  IDATA  changes  each  time  the  computer  passes  tnr  :uc 
this  part  of  the  program,  but  the  contents  of  array  IDATAT  is 
retained  and  added  to  each  time  the  computer  passes  through  tm 
part  of  the  program. 


IF  < RUN. EQ.  2 )  GO  TO  64® 

IF  ( RUN. EG. 3 )  GO  TO  665 
IF  ( RUN. EQ. 4 )  00  TO  69® 

IF  (RUN.EQ.5)  GO  TO  72® 

64®  CONTINUE 
REWIND  9 
REWIND  5 

READ (5) (IDATA <L>,L=1,N) 

READ ( 9 ) < I DAT AT ( LL ) , LL= 1 , NN ) 

DO  655  XX=1 , ELEM1 
DO  66®  YY=1 , 36 
ZZ=(XX-1 )*36 
DD=<  (PASS(XX)-l  )*13)+YY 
FF=YY+ZZ 

I DATAT ( FF ) = I DATAT ( FF ) + I DATA ( DD ) 
66®  CONTINUE 
655  CONTINUE 

WRITE(9) < I DAT AT ( LL ) , LL~ 1 , NN ) 

GO  TO  649 
665  CONTINUE 
REWIND  9 
REWIND  6 
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ije 


o  —  KtD  .  FOR 


F-  AD  <  4  )  ■  I DATA  <  L  1 ,  N ' 

-:t  AD  ■  -  >  i  I  DAT  AT  •  _L  )  .  L'_  - 1 ,  AM  ;> 

Jj  i.7*  XX~  1 ,  ELEi'll 

DO  o3U  VV-1,36 

2I~<  XX-1 >*36 

DD= ( ( PASS ( X  X ) - 1 ) *  1 3 ) +YY 

FF-YY+ZZ 

I DATAT ( FF ) = I DATAT ( FF ) + I DATA < DD ) 

680  CONTINUE 
670  CONTINUE 

WRITE ( 9 ) ( I DATAT ( LL ) , LL=1 , NN) 

GO  TO  649 
690  CONTINUE 
REWIND  9 
REWIND  7 

READ ( 7 ) ( I DATA (L) , L~1 , N) 

READ (9 ) < I DATAT ( L ) , LL= 1 ,  NN ) 

DO  700  X  X  =  1 , ELEM 1 

DO  710  YY~ 1 , 36 

ZZ=(XX-1 ) *36 

DD- < < PASS <  X  X ) - 1 ) * 1 3 ) +YY 

FF=YY+ZZ 

I DATAT  <  FF ) « I DATAT  <  FF ) + 1 DATA ( DD ) 

710  CONTINUE 
700  CONTINUE 

WRITE! 9) < I DAT AT  <  LL ) , LL= 1 ,  NN ) 

GO  TO  649 
720  CONTINUE 
REWIND  9 
REWIND  3 

READ ( 3 ) ( I DAT A  <  L ) , L= 1 , N  > 

READ ( 9 ) < I DAT AT ( LL ) , LL= 1 , NN ) 

DO  730  XX~1 , ELEIil 
DO  740  YY=1 , 36 
ZZ-=(XX-1  )*36 
DD=( (PASS(XX)-l )*13)+YY 
FF=YY+ZZ 

I DATAT ( FF ) = I DATAT ( FF ) + I DATA ( DD ) 

740  CONTINUE 
730  CONTINUE 

WRITE(9 ) ( I DATAT (LL ) , LL~1 , NN ) 

649  CONTINUE 

650  CONTINUE 
C 

-  We  now  have  array  IDATAT,  having  as  its  individual  elements 

-  the  total  of  five  individual  IDATA  elements  from  each  of 

-  five  elements.  To  get  the  average  value  of  these,  we  must 

-  divide  the  total  we  have  by  5.  This  division  is  done  in  the 

-  program  segement  below.  This  average  value  is  then  stored 

-  in  the  array  IDATAT  for  further  processing. 

REWIND  9 
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non  o  o 


F'  -i  q  i 


RED 


FOR 


REA  J  <  °  •  (  [  DA  :  A  V  i  L  >  ,  1_~  1  ,  i'ji'i  j 
DO  750  D I  --'—l « r:M 
I OATAT <DI7)~(  I D A ” A T  (  £'  1  V )  )  / 5 .  «i» 

750  CONTINUE 

WRITE (9) ( I DAT AT ( L ) , L= 1 , NN ) 

- —  The  steps  below  compute  airfoil  rotation  rate  (deg/sec), 

-  non-dimensional  rotation  rate,  Reynolds  number,  tunnel  "Q" 

-  and  volts  per  degree  for  the  run. 

REWIND  9 

READ ( 9 ) < IDATAT  <  L ) , L= 1 , NN ) 

DT I M=  < I DATAT  <  45 1 )-IDATAT(91 ) >*< 0.0010046) 

DPGSV - ((IDA  TAT ( 452 ) - 1 DATAT  <  92 ) > / 4096 . 0 ) *  1 0 . 0 
DPOSD~DPOSV  /  VF‘D 
RQTRAT  -=DPOSD/DT  I M 

NDRATE= ( R0TRA1 *.01 7453* 1 . 0 1 6 ) / ( 2 . 0*TUNVEL > 

RHO- ( BARQM*70 . 45 ) / (17 i 6 . 0* ( T60+TEMP ) ) 

MU~(2.270*(  10.0**(-3.0>  )*<  ( 460. 0+TEMP ) **1 . 5 )  )  /  (  460. 0+TEiv'F  + 1°8 ,  ; 
RE“  ( RHO  &TUNVEL*  1.01 6  )  /MU 
TUNQ= ( 0 . 5 ) *RHG* ( TUNVEL*  *2 > 

-  The  following  writes  pertinent  information  to  disk  file 

-  REDUDATADAT  as  a  heading. 

CALL  OPEN  <10,  -•  REDUDATADAT  '  ,  2 ) 

WRITE  (10,800) 

300  FORMAT  ( •’  DAY-1  ,  10X,  ''MONTH-'  ,9X,  'YEAR-'  ,9X,  •'  TIME'  ) 

WRITE  (10,810) DAY , MONTH , YEAR , HOUR 
310  FORMAT  ( 13, 11X, 13, 11X, 13, 9X, 15,/) 

WRITE  (10,820) 

320  FORMAT  ('  TEMPERATURE  '  ,  1 4X ,  '  BAROMETER'1  ) 

WRITE  ( 10,83@)TEMP,BAROM 
330  FORMAT  ( 2X, F6. 1 , 13X , F7. 2, / ) 

WRITE  (10,340) 

340  FORMAT  (  '  MANOMETER  1 ',  22X ,  ' MANOMETER  2"'  ) 

WR ITE  (10, 845 ) MANOM 1 , MAN0M2 
345  FORMAT  ( 2X , F8. 4 , 25X , F3. 4, / ) 

WRITE  (10,850) 

350  FORMAT  ( ■'  TUNNEL  VELOCITY' , 22X ,' MOTOR  VOLTAGE') 

WRITE  <10, 355 ) TUNVEL , MOTVOL 
355  FORMAT  ( 4X , F7. 2, 31 X , F6. 2, / ) 

WRITE  (10,360) 

860  FORMAT  ( '  DEGREES/SECOND' , 15X ,' NON-DIMENSIONAL  RATE ' ) 

WRITE  ( 10,370)ROTRAT,NDRATE 
370  FORMAT ( 4X , F6. 2, 25X , F8. 7, / ) 

WRITE  (10,330) 

330  FORMAT  ( •*  REYNOLDS  NUMBER' , 25X, ' TUNNEL  "Q"') 

WRITE  <10, 390 )  RE ,  TUNC! 

390  FORMAT  <  4X , El  1 . 4 , 30X , F6. 3, / ) 

DO  895  HH=1,16 

WRITE  (10, 397 ) HH, AVSTAT < HH ) 
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RED 


.  f  Ur 


120 


p  U  A  T"  1  :Vi- > Z-Ru—  iNF'L-T  rcMiu'^'A  :  .V'-ii DU  ■*  ^  -■  — 

C  DAT  ..  !.IUE 

One  pass  through  she  DO  ICO  J- 1 ,  N . 13  loop  computes  on  a 

point  in  the  CN  (normal  force  coefficient)  versus  ALPHA  curve . 

REWIND  9 

READ ( 9 ) ( I DATAT ( L ) , L= 1 , NN  > 

DO  100  J= 1 , NN, 36 
REDAT ( 1 )  =  I DAT AT  <  J ) 

REDAT ( 2 )  =  1 DAT  AT ( J+ 1 ) 

REDAT ( 1 9 ) = T DATAT ( J+l 8 ) 

REDAT ( 20 ) = I DATAT ( J+ 1 9 ) 

The  loop  below  subtracts  the  average  zero  input  readings 
(AVSTAT)  from  each  appropriate  IDATAT  element. 


DO  120  1=1,1:. 

AA-i J +1 ) +1 

REDAT i I+U  >  =1  DATA  I  ( AA ) - ( ( AVSTAT (I) ) -204 8 , 0 ) 
REDAT  (  I  +20 )  ~  I  DA"' AT  <  AA+ 1 8 )  -  (  ( AVSTAT  (  I  )  ) -2048.0) 
CONTINUE 


-  Operations  in  the  following  loop  correct  for  the  finite 

-  time  between  samples  using  a  linear  interpolation.  Time 

-  between  passes  must  be  sufficiently  small  or  the  linear 

-  interpolation  will  be  invalid. 

DO  140  R= 1 , 1 3 

REDATC ( R ) “REDAT <  R+ 1 3 ) - ( ( REDAT ( R+ 1 8 ) -REDAT <  R ) ) *( <  R-l ) / IS. 0 ) ) 
140  CONTINUE 


-  The  following  loop  converts  digital  quantities  to  degrees 

-  (angle  of  attack)  and  psi  (sensed  differential  pressure). 

-  The  AOA  conversion  below  assumes  the  A/D  board  is  strapped 

-  for  the  0—10  volt  unipolar  input  range.  The  amp  on  the 

-  board  is  set  for  a  gain  of  1,  so  any  input  to  the  board 

-  greater  than  10  volts  will  saturate  the  A/D  conversion  system. 

AOA=(  (  (REDATC(2>/4096.0)*10.0>— P0)/VF'D 
DO  160  S=l,16 
BB=( J+S)+l 

-  The  PRESS  conversion  below  assumes  the  A/D  board  is  strapp 

-  for  the  <-5)-(+5>  volt  bipolar  input  range,  where  the  inpu 

-  (from  the  transducers)  is  first  amplified  through  an 

-  amplifier  of  gain  100.  So  any  input  greater  than  +/-S0  mi i li- 

- volts  i.u ill  saturate  the  A/D  conversion  system. 


PRESS  <  S )  =  ( ( REDATC ( S+2 ) -2048 . 0 ) / 2043 . 0 ) *50 , © / SENS ( S ) 
CP  <  S )  =  ( PRESS ( S )  +  ( MANOM 1 / 27 . 68 ) ) / ( MAN0M2/ 27 . 68 ) 
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O  CONT  ~i  S'*  I !J !;£ 

-  '"hi  next  loop  defines  the  pressure  c  i  s  sr  iou  1 1  •:  n  on  --o  tjrn.;- 

C - surface  of  the  airfoi  I,  lead  in  a  ecae  to  tra  i  !  inq  ?c:ce. 

C  -  Pressure  coefficient  is  assumed  to  be  zero  at  the  trailing  edc. 

C 

WRITE  ( 10, 170)REDAT( 19) 

170  FORMAT  (///,-' TIME  BACK  TO  WHICH  DATA  FOR  THIS  PASS  HAS  BEEN 
+  CORRECTED: ' , F6. 0, /// ) 

WRITE  (10,185) 

185  FORMAT  ('  UPPER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L.E.  TO  T.E.,  ARE  GIVEN  BELOW',/) 

DO  200  V=l,9 
CPU ( V ) =CP ( V  > 

200  CONTINUE 

CPU ( 1 0 ) “CPU ( 9 )  +  ( 0 . 34 1 * ( CPU <  9 ) -CPU ( 8 ) ) ) 

DO  195  V=l,10 

WRITE  <1 0 , 1 90 ) V , CPU ( V ) 

1 90  FORMAT  ( •'  CPU  13,  FS.  4  > 

195  CONTINUE 
C 

-  The  next  loop  defines  the  pressure  distribution  on  the  lower 

-  surface  of  the  airfoil,  leading  edge  to  trailing  edge. 

-  Pressure  coefficient  is  assumed  to  be  zero  at  the  trailing  edge 

WRITE  (10,205) 

205  FORMAT  <//,•’  LOWER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L.E.  TO  T.E.,  ARE  GIVEN  BELOW’,/) 

CPL( 1 >=CP( 1 ) 

DO  220  W=2 , 8 

DD=18-W 

CPL ( W ) =CP ( DD ) 

220  CONTINUE 

CPL ( 9 ) =CPU ( 1@) 

DO  215  UNI, 9 
WRITE  ( 10,225)W,CPL(W) 

£25  FORMAT  ( ■’  CPL’ , 13, ' , F8. 4 ) 

215  CONTINUE 

—  The  following  loop  integrates  the  upper  pressure 
—  distribution  using  the  trapezoidal  rule. 

AREAUT=0.0 
DO  240  X=1 , 9 

LNGTHU=PORTU ( X+ 1 ) -RORTU ( X ) 

IF  ( (ABS(CPU(X+1>-CPU(X> ) >.GT. <ABS< <0.01 >*CPU<X) > ) )  GO  TO  245 
AREAU=  <0. 5)*< CPU<  X+l  )+CF’U<  X  >  >*LNGTHU 
245  IF  ( ( ABS ( CPU <  X  + 1 ) -CPU ( X ) ) ) . LE . <  ABS ( < 0 . 0 1 ) *CPU <  X ) ) ) )  GO  TO  235 
I NTU= ( PORTU  <  X ) -PORTU ( X+ 1 > ) *CPU ( X ) / ( CPU <  X  + 1 ) -CPU ( X ) > 

IF  ( I NTU . LT . LNGTHU )  GO  TO  260 
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V  —  RED  .FOR 

AREAU=  (  .  H  )  *  <  CPU  (  X  + 1 )  +CPIJ  (  X  )  )  *LNGTHU 
I™  ( I NTU ) .  GE .  (LiMGTHU)  >  CO  TO  235 

;  P  E  AU-  <  < .  5  >  *  I  NTU  --  r.FIJ  (  X  .  + 

<  <  .  3  >  <  (  Li'IQT HU-  i  N~ !j  )  <  OF  i_:  (  X  -1-  )  ) 

A  c  A  U  T  —  A  K  c.  A  U  T  ■*"  (  :X  E  A  U 

cOixt  i  nue 

o - The  following  loop  integrates  the  lower  pressure 

C  -  distribution  using  the  trapezoidal  rule. 

C 

AREALT =0 . 0 
DO  280  Y=1 , 8 

LNGTHL=PQRTL ( Y+ 1 ) -PORTL ( Y ) 

IF  (  ( ABS ( CPL <  Y+ 1  ) -CPL ( Y ) ) ) . GT . ( ABS ((0.01) *CPL ( Y ) ) ) )  GO  TO  270 
AREAL=  ( .  5 )  *  (*CPL  ( Y+ 1 )  +CPL  <  Y  )  )  *LNGTHL  " 

IF  ( < ABS ( CPL ( Y+ 1 ) -CPL ( Y ) ) ) . LE.  C ABS ( ( 0 . 0  1 ) *CPL ( Y ) ) ) )  GO  TO  275 
270  I NTL- ( PORTL ( Y ) -PORTL ( Y+ 1 ) ) *  CPL <  Y ) / ( CPL ( Y+ 1 ) -CPL ( Y  > ) 

IF  ( ( I NTL ) . LT. ( LNGTHL ) )  GO  TO  300 
AREAL- ( . 5 ) * ( CPL ( Y+ 1 ) +CPL ( Y ) ) «  LNGTHL 
IF  (< I NTL) .GE. (LNGTHL) )  GO  TO  275 
300  AREAL- ( ( . 5 ) *  I NTL*CPL ( Y ) )  + 

+  ( ( . 5  >  * ( LNGTHL- 1 NTL ) *CPL ( Y+ 1 ) ) 

275  AREALT =AREALT +AREAL 
280  CONTINUE 

NQRMCQ-AREALT-AREAUT 

C 

WRITE  (10, 900 ) AOA 

900  FORMAT  {>  ANGLE  OF  ATTACK--'  ,  F7.  3,  '  DEGREES') 

WRITE  (10,910) NORMCO 

910  FORMAT  (  •'  NORMAL  FORCE  COEFFICIENT- •'  ,  FS.  5,  /  ) 

WRITE  (10,950) 

950  FORMAT  ( •'  a*******************************-.'*********'#*****-*-#!** 

+********  fr*******--  ///  ) 

100  CONTINUE 
STOP 
END 
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=ROGRAM  DEN 
IMPLICIT  INTEGER  <A-Z) 

I N  TECER  I  DAT  A  <  505©  ) ,  i  DA  fA  1  (  35  i  -  [  DAT  AD  %  2  90 >  ,  PA’r  3  <  55 
I  NTEGER  N ,  R ,  X ,  V ,  V ,  W .  S ,  I ,  J ,  L ,  A  A ,  PP .  QQ 

INTEGER  RR,  SS,  TT,  IJIJ,  W,  WW.  XX,  YY,  22,  RUN,  TRAP ,  PAiZ ,  DIG,  NVMEL 
I NTEGER  ELEM 1 , ELEM2 , DAY , MONTH , YEAR , HOUR , CHANG 1 , CHANG2 
I NTEGER  DD , EE , FF , HH , LL , NN , BEG , START , END , FWD , HLT, MONO 
REAL  PORTLK  10) ,PORTL< 10) ,CP< 16) , CPU ( 10) , CPL< 10) , SENS< 16 ) 
REAL  IDATAT( 1260) , REDATC ( 36 ) , REDAT ( 36 ) , PRESS ( 16) 

REAL  P90 , P0 , TEMP , BAROM , MANOM 1 , MAN0M2 

REAL  TUNVEL, MOTVOL, A VST AT ( 16 ) , AREAUT, AREALT, RE, RHO, MU 
REAL  VPD, AOA, INTO, I NTL , NORMCG , TUNQ , LNGTHU , LNGTHL 
REAL  AREAU, AREAL, DTIM, DPOSD, DPGSV, ROTRAT, NDRATE 
REAL  AOASEP , RAWSEP , H I GH , LOW , TAKE 


- Load  transducer  sensitivities  ( mi  1 1 i vo 1 ts/os  i  ) 

DA  r A  SENS / 206 .5,174.2,174.9. 237 .5,207. 0 , 205 . 2 . 
►221.2, 172.3, 113. 1 , 1 13.3, 111.5, 138.7, 177.5,221 . 9 


-  Load  transducer  locations  on  upper  surface  (percent  chord* 

DATA  PORTU/0. 0, @. 0242,0. 0484 , 0. 0969, 0. 129, 0. 194, 0. 323 , 0. 60 
+0.888,1.000/ 

-  Load  transducer  locations  on  lower  surface  (percent  chord) 

DATA  PORTL/0. 0, 0.0161 , 0. 0319,0. 0484,0.0969,0. 194, 0. 323, 

+0. 686, 1 . 000/ 

WRITE  (1,5) 

5  FORMAT  (///,'  *****THE  DATA  FILES  TO  BE  REDUCED  MUST  BE  ON 
+  DISK  DRIVE  B  AND  MUST  BE  NAMED*****'’  ) 

WRITE  (1,6) 

6  FORMAT  ('  *************RAWDATA0. DAT.  RAWDATA1.DAT . 

+  RAWDATA5 .  DAT************  *•',///) 

-  Read  raw  data  from  RAWDATA0DAT  on  drive  B. 

CALL  OPEN  ( 3 ,  RAWDATA0DAT  ’  ,  2 ) 

READ  <3,1 0 ) DAY , MONTH , YEAR , HOUR 
10  FORMAT  (/, 13, 11X, 13, 11X, 13, 9X, 15) 

READ  ( 3 , 20 ) TEMP , BAROM 
20  FORMAT  < // , 2X , F6. 1 , 18X , F7. 2) 

READ  (3,30) MANOM 1 , MAN0M2 
30  FORMAT  <//,2X,F3.4,25X,F8.4> 

READ  (3, 40) TUNVEL, MOTVOL 
40  FORMAT  <//,4X,F7.2,31X,F6.2> 

READ  <3,50)P90,P0 
50  FORMAT  ( // , 5X , F7. 4, 23X , F7. 4 ) 

READ  (3,60) KOUNT , N 
60  FORMAT  <///,3X, I6,26X, 16) 

READ  <  3 , 70 ) AVSTAT ( 1 ) , AVSTAT ( 2 ) , A VST AT <  3 ) , AVSTAT ( 4 ) 

READ  <3,75  VAVSTAT <  5 ) , AVSTAT ( 6 ) , AVSTAT ( 7 ) , AVSTAT ( 8 ) 

READ  <  3, 75 ) AVSTAT  <  9 ) , AVSTAT ( 10) , AVSTAT (11), AVSTAT ( 12 ) 

READ  ( 3, 75) AVSTAT < 13) , AVSTAT ( 1 4 ) , AVSTAT ( 15 ) , AVSTAT ( 16) 

70  FORMAT  (////, F9. 3, 5X,  F'9.  3,  5X ,  F9.  3,3X,F9.3) 
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75  FORMAT  ■  F  9 . 3 , 5  X  ,  F  9 »  3 . 5  X  .  P  9 . 3 , 5  >  ,i-  .  3  ) 

VPD~  (  P 'V 0 **• P i  _> ;  /  ’”;0 . ft 

'0 - The  fol  i  owing  part  of  r  he  program  fir-zz  re  a  Is  the  or-jo  c  eg 

C - separation  AOA  ,  input  from  the  terminal.  I  he  corresponding 

C  -  raw  digital  separation  AOA  is  then  computed  and  used  as  follows 

C - LOW  is  computed  as  shown  be  1 ow,  then  used  in  the  program 

C  -  segement  beginning  with  EEG=2,  to  determine  the  first  AOA 

C  -  data  element  that  should  be  read  from  IDATA  in  RAWBATA1.DAT. 

C  -  This  AOA  data  element  is  designated  START,  and  the  AOA  data 

C  -  element  30  data  sets  later  is  designated  END.  With  the  first 

C  -  and  last  AOA  elements  known,  every  AOA  data  element  an  between 

C  -  read,  and  then  stored  in  the  IDATA1  array  for  matching  later  on 

C 

WRITE  < 1 , 85 > 

85  FORMAT  < '  ENTER  PREDICTED  SEPERAT  ION  AOA  (DEGREES)'  ,/> 

READ  ( 1 , 95 ) AOASEP 
95  FORMAT  (F5.2) 

RAWSEP-  (  (  ( AOAS£F’*VPD )  +P0 )  / 1 0 . 0 )  *4096 . 0 

IF  (<  MOTVOL .  GT .  5 .  00 ) .  AND .  <  MOT  VOL .  LE .  7 . 00 )  )  T  AKE~9<<> . 0 

IF  ( ( MOTVOL . GT . 7 . 00 ) . AND . ( MOTVOL . LE .  9 . 00 ) ) TAKE- 1 40 . 0 

IF  ( ( MOTVOL . GT . 9 . 00 ) . AND . ( MOTVOL . LE . 11 . 00 ) ) T AKE--200 . 0 

IF  <  <  MOTVOL .  GT .  1 1 . 00 ) .  AND .  ( MOTVOL .  LE .  1 3 .  00  )  >  TAKE-* 250 .  0 

LQW=RAWSEP~T  AKE 

CALL  0PEN<4, 'RAWDATA1DAT' ,2) 

READ ( 4 ) ( I DAT A ( L ) , L= 1 , N ) 

C 

BEG«2 

400  CONTINUE 

I F ( I DATA ( BEG ) . LE . LOW ) GO  TO  420 
START=BEG 
GO  TO  410 
420  BEG=BEG+13 
GO  TO  400 
410  CONTINUE 

END=3T  ART +540 
C 

ELEM1=0 

DO  500  PP"=START ,  END,  18 
ELEM 1 =ELEM 1 + 1 
I DATA 1 (ELEM1 >=IDATA(PP> 

500  CONTINUE 
C 

C  -  The  DO  650  loop, below  begins  by  branching  out  to  one  of  4 

C  -  program  locations  depending  on  the  number  of  times  the 

C  -  computer  has  been  through  the  650  loop.  The  program  segementc 

C  -  beginning  at  510,525,535  and  545  all  perform  the  same 

C  -  operation,  but  do  so  on  different  data  files,  in  turn. 

C  -  The  operation  performed  in  all  4  segements  is  as  follows. 

C  -  The  loop  contained  in  each  segement  reads  the  raw  digital 

C  -  angle  of  attack  from  every  data  set  in  the  array  IDATA. 

C  -  These  data  sets  were  previously  stored  in  the  second  through 
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FOR 


510 


52<d 


530 


540 

545 


550 


”  T  t  n  d  i  t  1  r  1  1  if  s  ,  .~’1  A I  *D  A  j  0  iL.  0;  1 1  t  r  r  >  tj  q  h  f '  A  -!lI  n  *  h5  '  .<  '  t  "i  i  "  i  -  ?_ 

oat  her  in  a  proqrarn  BIO.  Afti-r  z'-.nz~-  ancle:  o~  attic'-  =.* 
from  the  10ATA  array  i  •’<  the  proper  oata  rile,  the  i;  are  sc-:  r  ;q 
in  the  array  IDA7h2.  Note  that  the  contents  of  the  /:?;“>  7"^  s.-:-- 
are  replaced  each  time  the  DO  650  loop  it  travelled. 

DO  650  RUN=2, 5 
ELEM2=0 

IF  (RUN.EQ.2)  GO  TO  510 

IF  (RUN.EQ.3)  GO  TO  525 

IF  (RUN. EQ. 4 )  GO  TO  535 

IF  (RUN.EQ.5)  GO  TO  545 

CONTINUE 

CALL  OPEN  ( 5 .  RAWDATA2GAT  1 , 2  > 

READ ( 5 ) ( I DAT A  <  L ) , L= 1 ,  N ) 

DO  520  QQ-2 , N , 1 8 

ELEM2=ELEM2+1 

I DAT A2 ( ELEM2 ) ~ I DATA  <  QQ ) 

CONT I NUE 
GO  TO  550 
CONTINUE 

CALL  OPEN (6, 1 RAWDATA3DAT J  ,  2  > 

READ ( 6 ) ( I DATA (L ) , L=1 , N) 

DO  530  QQ=2,N,18 

ELEM2=ELEM2+1 

I DATA2 ( ELEM2 ) = I DATA ( QQ ) 

CONTINUE 
GO  TO  550 
CONTINUE 

CALL  OPEN ( 7 ,  RAWDATA4DAT '  ,2) 

READ ( 7 ) < I DAT A ( L ) , L= 1 , N ) 

DO  540  QQ=2, N, IS 

ELEM2=ELEM2+1 

I DATA2 ( ELEM2 ) = I DATA  <  QQ ) 

CONTINUE 
GO  TO  550 
CONTINUE 

CALL  OPEN < 8 ,  RAWDATA5DAT  ‘  ,2) 

READ ( 8 ) < I DATA ( L ) , L= 1 , N ) 

DO  550  QQ— 2, N, 18 

ELEM2=ELEM2+1 

I DATA2 ( ELEM2 ) = I DAT A ( QQ ) 

CONTINUE 


-  We  now  have  an  array,  IDATA1,  containing  every  AOA  element 

-  from  within  a  specified  interval  in  the  array  IDATA, 

-  as  read  from  disk  file  RAWDATA1.DAT.  We  also  have  an 

-  array,  IDATA2,  containing  every  angle  of  attack  element 

-  from  one  of  the  four  remaining  disk  files.  The  program  segemen 

-  below  compares  each  element  in  IDATA1  with  every  element  in 

-  1DATA2,  in  turn.  When  the  closest  match  it  found,  this  infor- 

-  rnation  is  stored  in  the  PASS  array  for  use  in  the  next  pro- 
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■- - ora r .i  s <t  q  =•  r»i <=» r i . 

C 

RR=-£> 

DO  560  SS=1,ELEM1 
CHANG 2® 4 09 & 

TRAP*® 

DO  57®  TT=1  f  ELEM'2 
RRaRR+1 

CHANG1=I ABS ( I DATA 1 ( SS ) -IDATA2<  TT ) ) 

IF  < CHANG 1.GE.CHANG2)  GO  TO  580 
CHANG2=CHANG1 
GO  TO  570 
580  TRAP=TRAP+1 
IF  ( TRAP. NE. 1 )  GO  TO  570 
PAS'S  ( SS )  =RR- 1 
570  CONTINUE 
560  CONTINUE 

C  -  The  next  part,  up  to  and  including  line  610,  reads  every  data 

C - set  from  I DA  FA,  in  disk  file  RAWDATA1.DAT,  that  corresponds  t 

C  -  an  AOA  element  contained  in  array  IDATA1.  That  is,  every  data 

C  -  set  within  the  previously  defined  interval  is  read  from  the 

C - 1DATA  array  in  disk  file  RAWDATA1.DAT.  These  data  sets  are  tl< 

C  -  repacked  sequentially  in  the  array  IDATAT. 

C 

IF  ( RUN. NE. 2 )  GO  TO  630 
REWIND  4 

READ ( 4 ) < IDATA(L) ,L=1,N> 

CALL  OPEN  ( 9 ,  AVRGDA  T ADAT  %  2 ) 

PAZZ~0 

FWD=START-1 

HLT=END-1 

DO  610  UU=FWD , HLT , 36 
PAZZ=PAZZ+1 
DO  620  VV~ 1 , 36 
DD=(UU+VV)-1 
WW=VV+ ( 36* ( PAZZ-1 ) ) 

IDATAT (WW)=IDATA(DD) 

620  CONTINUE 
610  CONTINUE 

NN=PAZZ*36 

WRITE (9) ( I DAT AT ( LL ) , LL= 1 , NN ) 

630  CONTINUE 


-  The  following  program  segements  beginning  on  lines  640,665, 

-  690,  and  720  perform  the  same  operation,  but  does  so  on 

-  different  disk  files,  in  turn.  The  operation  performed  is  as 

-  follows.  The  information  previously  obtained  and  stored  in 

-  the  array  PASS  is  now  used  to  determine  which  data  sets  will 

- be  read  (from  whichever  disk  file  of  the  remaining  four-  that 

-  is  then  being  worked  on),  then  added,  element  by  element,  to 

-  the  element  totals  in  the  array  IDATAT.  Note  that  the  content 
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C  -  of  tr*  arra^  IDATA  changes  eaO  tins  the  cor.  cut  er  * 

C  -  trie  pa-t  of  the  croc  ram,  rut  the  ;c  n terra  c'  a.---n  i  FD- 

C - ratain-i'j  ano  added  to  each  t i rue  the  comouta1--  oa:r.  v  ;  nv 

0  -  part  or  the  program. 

C 

IF  (RUN. EQ. 2)  GO  TO  640 
IF  (RUN.EQ.3)  GO  TO  665 
IF  ( RUN. EQ. 4 )  GO  TO  690 
IF  ( RUN. EQ. 5 )  GO  TO  720 
640  CONTINUE 
REWIND  9 
REWIND  5 

READ  <  5 ) < IDATA(L) , L=1 , N) 

READ ( 9 ) ( I DATAT  <  LL ) , LL- 1 , NN  > 

DO  655  XX-1 , ELEM1 

DO  660  YY  = 1 , 36 

ZZ=<  XX-1 ) *36 

DO- ( <  PASS <  X  X ) - 1 ) * 1 8  >  +YY 

f F-YY+ZZ 

I DATAT ( FF )  =  I DATAT ( F  F )  + 1 DATA ( DD ) 

660  CONTINUE 
655  CONTINUE 

WRITE(9) ( I DAT  AT ( LL ) , LL® 1 ,  NN ) 

GO  TO  649 
665  CONTINUE 
REWIND  9 
REWIND  6 

READ (6) < I DATA ( L ) , L-l , N ) 

READ(9) ( I DAT AT ( LL ) , LL- 1 ,  NN ) 

DO  670  XX=1 , ELEM1 
DO  680  YY— 1 , 36 
ZZ-(XX-1 )*36 
OD=( (PASS(XX)-l )*13>+YY 
FF-YY+ZZ 

I DATAT ( FF ) - 1 DATAT  <  FF  >  + 1 DA T A ( DD ) 

630  CONTINUE 
670  CONTINUE 

WRITE ( 9) ( I DATAT ( LL ) , LL=1 , NN) 

GO  TO  649 
690  CONTINUE 
REWIND  9 
REWIND  7 

READ(7) ( IDATA(L) ,L=1 ,N) 

READ (9X1  DATAT (L) ,  LL- 1 ,  NN ) 

DO  700  XX-1 , ELEM1 

DO  710  YY— 1 1 36 

ZZ=(XX-1 )*36 

DD® ( ( PASS ( X  X ) - 1 ) * 1 8 ) + YY 

FF-YY+ZZ 

I DATAT ( FF )  =  I DATAT ( FF  >  + 1 DATA  <  DD ) 

710  CONTINUE 
700  CONTINUE 
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toft  I  Ts  <  9  ;  I  OR  1  AT  \  LL )  ,  LL=  1 ,  NN  > 

GO  TO  649 
72'o  CONTINUE 
REWIND  9 
REWIND  8 

READ  <  8  > ( IDATA(L) , L=1 ,N) 

READ (9) ( I DATAT <  LL ) , LL= 1 , NN  > 

DO  730  XX=1,ELEM1 

DO  740  YY=1 , 36 

22= ( XX-1 )*3 6 

DD= ( ( PASS ( X  X ) - 1 ) * 1 8 ) +YY  ' 

FF=YY+ZZ 

I DATAT ( FF ) = I DATAT ( FF  >  + 1 DATA ( DD  > 

740  CONTINUE 
730  CONTINUE 

WRITE  <  9) ( I DAT AT ( LL ) , LL= 1 ,  NN ) 

649  CONT I NUE 

650  CONTINUE 
C 

C  -  We  now  have  array  1 DATAT,  having  as  its  individual  element: 

C  -  the  total  of  five  individual  IDATA  elements  from  each  of 

C  -  five  elements.  To  get  the  average  value  of  these,  we  must 

C  -  divide  the  total  we  have  by  5.  This  division  is  done  in  the 

C  -  program  segement  below.  This  average  value  is  then  stored 

C  -  in  the  array  IDATAT  for  further  processing. 

C 

REWIND  9 

READ ( 9 ) ( IDATAT ( L ) , L=i , NN ) 

DO  750  D1V=1,NN 
IDATAT < DIV )  =  ( IDATAT(DIV) )/5. 0 
750  CONTINUE 

WRITEI9) ( I DAT AT ( L ) , L= 1 , NN  > 

C 

C  -  The  steps  below  compute  airfoil  rotation  rate  (deq/sec), 

C  -  non-dimensional  rotation  rate,  Reynolds  number,  tunnel  "0" 

C  -  and  volts  per  degree  for  the  run. 

C 

REWIND  9 

READ<9) ( IDATAT (L) , L=1 , NN) 

DTIM=< IDATAT (451 ) -IDATAT (91 ) >*(0.0010046) 

DPQSV=< ( IDATAT (452) -IDATAT (92) > M096. 0 >*10. 0 

DPOSD=DPOSV/VPD 

R0TRAT=DP03D/DT I M 

NDRATE= ( ROTRAT* . 0 17453* 1 . 0 1 6  > / ( 2 . 0*TUNVEL  > 

RHO= ( BAROM*  70 . 45 ) / ( 1 7 1 6 . 0* ( 460+TEMP  > ) 

MU= (2. 270* ( 10.0**<-8.0>  >*( ( 460 . 0+TEMP )**1.5) ) / (460.0+TEMP+193. 
RE= ( RHO*TUNVEL* 1 . 0 1 6 ) / MU 
TUNQ= (0.5) *RHO* ( TUNVEL**2  > 

C 

C  -  The  following  writes  pertinent  information  to  disk  file 

C  -  REDUDATADAT  as  a  heading. 

C  102 
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CALL  OPEN i  0 ,  •'  PEDUDATADA  i  ,  L 
WRIYE  ( 10,  300  '• 

S00  FORMA  r  ‘  DAY  1  ,  10X, -1  MONTH YEAK1  ,  vX  1  IM£-'  • 

(*i  RE  i  E  (1 0 ,  8 1  @  >  D  A  Y ,  MONTH,  YEAR,  H  0  U  R 
810  FORMAT  (13, 1 1 X , 13, 1 1 X , I 3, 9X , 15,/) 

WRITE  (10,320) 

320  FORMAT  ('  TEMPERATURE' , 1 4X BAROMETER'  ) 

WRITE  ( 10,330>TEMP,BAROM 
830  FORMAT  (2X,F6. 1 , 18X, F7. 2, / ) 

WRITE  (10,340) 

840  FORMAT  ('  MANOMETER  1 1 , 22X , ' MANOMETER  2 ') 

WRITE  ( 10,345)MANOM1,MANOM2 
845  FORMAT  (2X , F3. 4, 25X , F8. 4, / ) 

WRITE  (10,850) 

350  FORMAT  ( •'  TUNNEL  VELOCITY1  , 22X ,  1  MOTOR  VOLTAGE1) 

WR I TE  (10, 355 ) TUNVEL . MOT VOL 
855  FORMAT  ( 4 X , F7 . 2 , 3 1 X , F6 . 2 , / ) 

WR I TE  ( 10, 860 ) 

860  FORMAT  (  -  DEGREES/ SECOND .  15X,  •' NON-DIMENSIONAL  RATE') 

WRITE  <1 0 . 870 ) ROTRAT , NDRATE 
370  FORMAT ( 4  X , F6 . 2 , 25X , F3 . 7 , / ) 

WRITE  (10,380) 

830  FORMAT  ( •’  REYNOLDS  NUMBER 1 , 25X ,  1  TUNNEL  "O'"  ) 

WRITE  ( 10, 890) RE, TUNG 
390  FORMAT  <4X , El 1 . 4, 30X , F6. 3, / ) 

DO  395  HH-1,16 

WRITE  (10, 897 ) HH , AVST AT ( HH ) 

397  FORMAT  ( '  AVERAGE  ZERO- INPUT  READING,  TRANSDUCER ' , 13, '  ~',F  . 
395  CONTINUE 
C 

C  -  One  pass  through  the  DO  100  J= 1 , NONO , 1 8  loop  computes  one 

C  -  point  in  the  CN  (normal  force  coefficient)  venous  ALPHA  cur 

C 

REWIND  9 

READ ( 9 ) ( I DAT AT ( L  > , L= 1 , NN ) 

NONO~NN- 1 3 

DO  100  >J=1 ,  NONO,  13 

REDAT ( 1 ) *  I DATAT ( J ) 

REDAT ( 2 )  =  I DAT AT ( U+ 1 ) 

REDAT  ( 19)*IDATAT  ( J+18) 

REDAT ( 20 ) » I DATAT  <  J+ 1 9  > 

- The  loop  below  subtracts  the  average  zero  input  readings 

-  (AVSTAT)  from  each  appropriate  IDATAT  element. 

DO  120  1*1,16 
AA*<  J+I ) +1 

REDAT ( 1+2) =1 DATAT (AA)-( ( AVSTAT ( I )) -2048.0) 

REDAT (I +20 >* I DATAT <  AA+18)-<  ( AVSTAT  ( I )  ) -2043.0) 

120  CONTINUE 

-  Operations  in  the  following  loop  correct  for  the  finite 
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-  -  tir \\z  beti'i-isn  samp  .  as  using  a  linear  1  n  Sirpc  ’a’ : or .  :  1  3 

-  bit'ijsen  passes  must  be  sufficiently  smal  i  •-•r-  ~h  a  linear 

o - inserpolati o n  w ill  be  i n v a  1  id. 

f: 

DO  140  R=l,13 

REDATC < R )  =REDAT ( R+18 ) - (  < REDAT ( R+ 1 8  > -REDAT ( R )  )*(  (R-D/13.0)  ) 

140  CONTINUE 

-  The  following  loop  converts  digital  quantities  to  degrees 

-  (angle  of  attack)  and  psi  (sensed  differential  pressure). 

-  The  AOA  conversion  below  assumes  the  A/D  board  is  strapped 

-  for  the  0-10  volt  unipolar  input  range.  The  amp  on  the 

-  board  is  set  for  a  gain  of  1,  so  any  input  to  the  board 

-  greater  than  1@  volts  will  saturate  the  A/D  conversion  system. 

AOA* < ( <  REDATC ( 2 ) / 4096 . 0 ) *  1 0 . 0 ) -P0 ) / VPD 
DO  160  S= 1 , 16 
BB~ ( J+S ) +1 

-  The  PRESS  conversion  below  assumes  the  A/D  board  is  strapned 

-  for  the  (-5)-(+5)  volt  bipolar  input  range,  where  the  input 

-  (from  the  transducers)  is  first  amplified  through  an 

- amplifier  of  gain  100.  So  any  input  greater  than  +  /~50  mil  I i  — 

-  volts  will  saturate  the  A/D  conversion  system. 

PRESS ( S )  =  ( ( REDATC ( S+2 ) -  2048 . 0 ) / 2048 . 0 ) *50 . 0 / SENS ( S ) 

CP  <  S )  =  ( PRESS ( S )  +  ( MANOM 1 / 27 . 63  > ) / ( MAN0M2 / 27 . 63 ) 

60  CONTINUE 

- The  next  loop  defines  the  pressure  distribution  on  the  upper 

-  surface  of  the  airfoil,  leading  edge  to  trailing  edge. 

-  Pressure  coefficient  is  assumed  to  be  zero  at  the  trailing  edge 

WRITE  ( 10, 1 70) REDAT ( 19) 

170  FORMAT  (//,' TIME  BACK  TO  WHICH  DATA  FOR  THIS  PASS  HAS  PEEN 
+  CORRECTED:  •' ,F 6.0, /) 

WRITE  (10,185) 

135  FORMAT  ( '  UPPER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L.E.  TO  T.E.,  ARE  GIVEN  BELOW1,/) 

DO  200  V= 1 , 9 
CPU(V)=CP(V) 

200  CONTINUE 

CPU( 10)=CPU(9>+(0.341*(CPU(?)-CPU<8) ) ) 

DO  195  V* 1,10 
WRITE  (10, 190)V,CPU(V) 

190  FORMAT  ('  CPU' , 13, ,F3.4> 

195  CONTINUE 

-  The  next  loop  defines  the  pressure  distribution  on  the  lower 

-  surface  of  the  airfoil,  leading  edge  to  trailing  edge. 

-  Pressure  coefficient  is  assumed  to  be  zero  at  the  trailing  edge. 
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W'K  l  TE  »  1 ,  205 ) 

..05  FORMAT  ( /  / ,  '  LONER  SURFACE  PRESSURE  COEFFICIENTS, 

+  L .  b. .  TO  T.L.,  hRE  GIVEN  BELOW',/) 

CF'L<  1  >=CP<  1  ) 

DO  220  W=2 , S 

DD=18-W 

CPL ( W ) =CP ( DD ) 

220  CONTINUE 

CPL ( 9 ) =CPU (1@) 

DO  215  W=1 , 9 

WRITE  ( 10,225)W,CPL(W) 

225  FORMAT  ( •*  CPL  * ,  13,  '=•’  ,F3.4) 

215  CONTINUE 
C 

C  -  The  following  loop  integrates  the  upper  pressure 

C  -  distribution  using  the  trapezoidal  rule. 

C 

AREAUT *0 .  <:> 

DO  240  X  =  1 , 9 

LNG THU=PORTU  <  X+ 1  )  -PORTU  (  X  > 

IF  ( <ABS( CPU (X+l >-CPU(X> )) .GT. (ABS < (0.01 > *CPU(X) >) )  GO  TO  215 
AREAU* ( 0 . 5 ) * ( CPU ( X+ 1 ) +CPU ( X ) ) *LNGTHU 
245  IF  ( (ABS(CPU(X+1 )-CPU(X> ) >.LE. <ABS( (0.01 >*CPU(X) > ) >  GO  TO  235 
INTU--  ( PORTU  (  X  )  -PORTU  (  X  + 1 )  )  *CPU ( X  )  /  ( CPU  ( X+ 1 )  -CPU  (  X  )  ) 

IF  (INTU.LT.LNGTHIJ)  GO  TO  260 
AREAU=  < . 5 ) * ( CPU ( X+l > +CPU ( X ) ) *LNGTHU 
IF  ( ( INTU ) . GE. ( LNGTHU ) )  GO  TO  235 
260  AREAU= ( ( . 5 ) * I NTU*CPU ( X ) ) + 

+  ( ( . 5 ) * ( LNGTHU- I NTU )  •*CPU  (  X  + 1 ) ) 

235  AREAUT  *AREAUT +AREAU 

240  CONTINUE 


The  following  loop  integrates  the  lower  pressure 
distribution  using  the  trapezoidal  rule. 


AREALT=0. 0 
DO  230  Y*1 , 8 

LNGTHL=PORTL <  Y+ 1  > -PORTL ( Y ) 

IF  ( ( AB3 ( CPL ( Y+ 1 ) -CPL ( Y ) ) ) . GT . ( ABS ((0.01) *CPL ( Y ) > ) )  GOTO  270 
AREAL* < . 5 ) * ( CPL ( Y+ 1 ) +CPL ( Y ) ) *LNGTHL 

IF  ( ( ABS ( CPL ( Y+ 1 ) -CPL  < Y ) ) ) . LE . ( ABS ( ( 0 . 0 1 ) *CPL ( Y ) ) ) )  GO  TO  275 
270  I NTL* ( PORTL ( Y ) -PORTL ( Y+ 1 > ) *CPL ( Y ) / ( CPL( Y+l ) -CPL ( Y ) ) 

IF  ( (INTL).LT. (LNGTHL) )  GO  TO  300 
AREAL* ( . 5 ) * ( CPL ( Y+ 1 ) +CPL <  Y ) > *LNGTHL 
IF  ( (INTL).GE. (LNGTHL) )  GO  TO  275 
300  AREAL* ( (.5)*INTL*CPL(Y) )+ 

+  ( (.5>*(LNGTHL-INTL)*CPL(Y+1> ) 

275  AREALT =AREALT +AREAL 

230  CONTINUE 

NORMCQ=AREALT-AREAUT 
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'  FORMAT  f  ANGLE  jF  AT  7AC<"-‘  ,  FT  ,5,  '  LEGh:E  ES  '  ) 

W-  i~t  i  >■:  ,  v  1  <o  >  (MOrtHCU 

9 1  t>  FORMAT  NORMAL  FORCE  COEFFICIENT ,  FS.  5,  /  ) 
WRITE  <10,950) 

950  FORMAT  (  1  ********************************** 

+  *******•»■»*•»•**■»■»*•’  ,  ///  ) 

100  CONTINUE 
STOP 
END 
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P R 0 GRAM  f  E ST 

P OR  CAL i SPAT iuN  OP  T HE  HIM  12  IN  3 i PuLAK .  10 

CHANNEL  USED  IS  ALWAYS  CHANNEL  l 

LINK:  I  ESI ...  ALIO,  FOKLIB/S,  TEST/N/E 

REAL  VALTOT 

IMPLICIT  INTEGER  (A-Z) 

WRITE  (1,20) 

20  FORMAT  ENTER  CHANNEL  YOU  WISH  TO  SEE  (0-1 

READ  (1,30) CHAN 
30  FORMAT  (14) 

10  CONTINUE 

VALTQT=0. 0 

DO  100  1=1,100 

CALL  AD (VALUE, CHAN, 80) 

VALTOT=VALTOT+( VALUE/ 100. 0) 

100  CONTINUE 

WRITE  (1,40)  CHAN-,  VALTOT 

40  FORMAT  (•-►AVERAGE  VALUE,  CHANNEL  ■'  ,  14,  •'  =  •'  ,  F 

GO  TO  10 
STOP 
END 
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1  — 

AD  1 0 

ivIAC 

.  --i'j 
d  T  r. \ 

HU 

H  /  U  C 

Eft-, 

'ICE  Ra 

time 

r  OR  f  P 

AM 

CALL ALL 

e 

CALL 

AD ( VALUE , C 

HAN ,  f.ASE ) 

GET  ONE 

SAMPLE 

FROM  THE  CHAN' TH  CHANNEL 

ON  THE  A/D  BOARD  WITH  BASE  ADDRESS  •’  BASE •' 

AD: 

LD 

( VALUE ) 

,  HL 

LD 

(CHAN), 

DE 

LD 

(BASE) , 

BC 

EX 

DE,  HL 

; HL->CHAN 

LD 

A, (HL) 

; GET  CHAN  NO. 

LD 

HL, (BASE) 

LD 

C,  (HL) 

; GET  BASE  I/O  ADDRESS  TO  C  REG  FOR 

OUT 

(C),  A 

; MODE  0  TO  CHAN  MO. 

; USES  BASE  ADDRESS  IN  C  R03 

INC 

c 

; POINT  TO  START  CONVERSION  PORT 

LD 

A,  O 

OUT 

(C)  ,  A 

; START  CONVERSION 

DEC 

c 

; POINT  TO  BASE  REGISTER 

NRDY: 

IN 

A,  (C) 

; GET  STATUS 

AND 

08®H 

; BIT  7  IS  STATUS,  =1  IS  BUSY 

JR 

N2, NRDY 

; NOT  ALL  0‘S  =>  BUSY 

INC 

C 

; POINT  TO  BASE  ADD+1 

INC 

C 

; POINT  TO  DRL 

IN 

A,  (C) 

; LOU  BYTE  OF  VALUE 

LD 

E,A 

INC 

C 

; POINT  TO  DRH 

IN 

A,  (C) 

; HIGH  BYTE  OF  VALUE 

AND 

0FH 

; MASK  OUT  HIGH  NIBBLE 

LD 

D,  A 

: DE=VALUE 

LD 

HL, (VALUE)  :HL-> WHERE  TO  PUT  VALUE 

LD 

( HL ) , E 

; PUT  LOW  BYTE  OF  VALUE 

INC 

HL 

LD 

(HL) ,D 

; THAT  GIVES  THE  CALLER  THE  VALUE 

? 

RET 

VALUE: 

DW 

0 

; STORAGE  FOR  ADDRESS  OF  VALUE 

CHAN: 

DU 

0 

; STORAGE  FOR  ADDRESS  OF  CHANNEL  NO 

BASE: 

DU 

0 

; STORAGE  FOR  ADDRESS  OF  BASE  ADDRES 

; 

.  Z80 

ENTRY  DA 

• 

* 

CALL 

DA (VAL, CHAN, BASE) 

• 

♦ 

CHAN 

IS 

0 

• 

♦ 

BASE 

IS 

72 (BASE  10) 

DA: 

LD 

A, (DE) 

; GET  CHAN 

ADD 

A,  A 

; DOUBLE  IT 

INC 

A 

; ADD  ONE 
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AD  10 


MAC 


PUSH 

si‘|  ICbJ 

HL 

POP 

BC 

HL 

LO 

C, ( HL ) 

ADD 

A,  C 

LD 

C,  A 

POP 

HL 

LD 

A, <  HL ) 

OUT 

<  C ) ,  A 

DEC 

C 

INC 

HL 

LD 

A,  (HL) 

OUT 

<C),A 

RET 

END 

;  SAVE  VAL 
‘  J,u“  ■  BASE 

•  C -LOW  t  Y T  t  OF  Bhs£ 

C-LOW  BYTE  VALUE  nF  part 

get  VAL 

GET  LOW  BYTE 

PUT  LOU  BYTE 

C=HIGH  BYTE  PORT 

HL=  >H I  BYTE 

GET  HI  BYTE 

PUT  HI  BYTE 


I 


Rage 


1 


ONE 


—  STCLK 


ENTRY  TCL-. 

L  ■  ;  LO  A.vl  7H  j  O' .-in.  i.“iE_  1  CTri.  1--L  “T.f".Z  . 
COT  >.  079H  )  ,  A 

L  Ci  A ,  o vAH  :  T 1  ME  Co  N3T AnT 
OUT  ( 079H ) , A 

LD  A, 057H  ; CHANNEL  2  CNTRL  WD=CTR 

OUT  ( 07AH ) , A 

LD  A, 0FFH  ; TIME  C0NSTANT=256(EASE  10) 
OUT  ( 07AH) , A 

LD  A, 057H  ; CHANNEL  3  CNTRL  WD=CTR 

OUT  <07EH),A 

LD  A, 0FFH  ;  TIME  CONSTANT 
OUT  ( 07BH ) , A 

JP  0  ; SYSTEM  REENTRY  POINT 

END  STCLK 
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FORCE  COEFFICIENT 


0.(H  '."t'-.T  ■  l-'-i  t  v-t.  f-  i--t  ■,  --.V . i—  f-:r- 

o  5  10  15  20  25  30  35 

ANGLE  OF  ATTACK  (DEGREES) 


Figure  21.  Force  Coefficient  vs.  Angle  of  Attack 

for  <*=32.8  Deg/Sec  and  \^a26.7  Ft/Sec 
(  aNDa0.0il  ) 
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10  15  20  25 

ANGLE  OP  ATTACK  (DEGREES) 


Figure  23.  Force  Coefficient  vs.  Angle  of  Attack 

for  6=74.5  Deg/Sec  and  ^=26.7  Ft/Sec 
(  ®nd*0.025  ) 


ANGLE  OF  ATTACK  (DEGREES) 


Figure  24.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =97.6  Deg/Sec  and  Vms26  .7  Ft/Sec 
(  ®nd“0*032  ) 


-5  0  5  lo  15  20 

ANGLE  OF  ATTACK  (DEGREES) 


Figure  25.  Force  Coefficient  vs.  Angle  of  Attack 
for  a  =0.0  Deg/Sec  and  X»=2 6.7  Ft/Sec 
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FORCE  COEFFICIENT 


15  20  25 

ATTACK  (DEGREES) 


Figure  26.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =33.8  De&/Sec  and  Yo=30.1  Ft/Sec 
(  ^ND=®  * oio  ) 


FORCE  COEFFICIENT 


Figure  27 


Force  Coefficient  vs.  Angle  of  Attack 
for  a  =54, 8  Deg/Sec  and  ^  =  30,1  Ft/Sec 
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1. 


FORCE  COEFFICIENT 


ANGLE  OF  ATTACK  (DEGREES) 


Figure  28.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =75 .6  Deg/Sec  and  Y=30.i  Ft/Sec 

(  <*nd=0.022  ) 


120 


FORCE  COEFFICIENT 


Figure  29  .  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =95.1  Deg/Sec  and  \  =  30. 1  Ft/Sec 
(  «ND=0.029  ) 
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Figure  31.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =32  .7  De^Sec  and  ^=39, 9  Ft/Sec 
(  ®j$ds0.008  ) 
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ANGLE  OF  ATTACK  (DEGREES) 


Figure  32.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =  55.8  Beg/Sec  and  Vro=39.9  Ft/Sec 

(  ) 


FORCE  COEFFICIENT 


Figure  33.  Force  Coefficient  vs.  Angle  of  Attack 

for  <*=75. 2  Deg/Sec  and  Tfc=39.9  Ft/Sec 
(  aND=0. 017  ) 
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FORCE  COEFFICIENT 


9 


Figure  34.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =97.1  Deg/Sec  and  V=39.9  Ft/Sec 
(  ®pju=0.022  ) 
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FORCE  COEFFICIENT 


Figure  35.  Force  Coefficient  vs.  Angle  of  Attack 
for  Q= 0.0  Deg/Sec  and  X»=39 *9  Ft/Sec 
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ANGLE  OF  ATTACK  (DEGREES) 


Figure  36.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =  32.  7  Deg/Sec  and  ^=47. 8  Ft/Sec 
(  *  006  ) 


FORCE  COEFFICIENT 


Figure  37.  Force  Coefficient  vs.  Angle  of  Attack 

for  fl*53.4  Deg/Sec  and  ^=47.8  Ft/Sec 

(  <*nd»0.010  ) 
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\ 


i 


ANGLE  OF  ATTACK  (DEGREES) 


Figure  38.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =74. 8  Deg/Sec  and  V  =47.  8  Ft/Sec 
(  «NDs0.014  )  °° 


Figure  39.  Force  Coefficient  vs.  Angle  of  Attack 

for  a  =97.1  Deg/Sec  and  Ft/Sec 

(  aNDa°.0l8  ) 


FORCE  COEFFICIENT 


Figure  40.  Force  Coefficient  vs.  Angle  of  Attack 
for  a  *0.0  Deg/Sec  and  ^»=47.8  Ft/Sec 
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